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ABSTRACT 


This  thesis  reviews  Runge-Kutta-Nystrom  (RKN)  methods 
for  the  solution  of  second  order  ODE's.  Presented  are 
new  derivations  of  the  constraints  for  methods  of  rank 
three  and  four  together  with  the  corresponding  truncation 
error  bounds. 

Two  procedures  are  employed  to  obtain  new  coefficients 
for  RKN  methods.  The  use  of  these  coefficients  is  evaluated 
with  regards  to  truncation  error  bounds,  accuracy  and  stab¬ 
ility.  Numerical  experiments  were  performed  on  an  IBM 
360/67  and  comparisons  are  made  between  the  results  of 
solving  a  second  order  equation  directly  by  use  of  RKN 
methods  and  by  an  RK  method  applied  to  the  equivalent 
system  of  first  order  equations.  It  is  shown  that  RKN 
methods  are  in  general  superior  to  RK  methods. 
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PREFACE 


This  thesis  reviews  Runge-Kutta  type  methods  for 
solution  of  second  order  ordinary  differential  equations. 
The  solution  of  such  an  ODE  by  a  simultaneous  system  of 
first  order  equations,  and  directly,  by  the  Runge-Kutta- 
Nystrom  (RKN)  method,  will  be  discussed. 

New  sets  of  coefficients  for  RKN  methods  will  be 
presented.  The  solution  of  a  second  order  equation  using 
these  methods  will  be  compared  with  the  solution  by  a 
system  of  equations  with  regards  to  accuracy,  stability 


and  truncation  error. 
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CHAPTER  I 


INTRODUCTION 


1. 1  The  Runge-Kutta  Method 

The  numerical  solution  of  differential  equations 
is  an  interesting  and  important  problem  in  the  area  of 
applied  mathematics.  In  particular  ordinary  differential 
equations  (ODE)  of  first,  second  and  higher  order  are 
often  encountered,  a  large  class  of  which,  though  linear 
and  of  a  simple  form,  cannot  be  solved  analytically. 

This  has  prompted  the  development  of  several  numerical 
methods  of  which  the  basic  idea  for  one  was  first  pro¬ 
posed  by  Runge  [20],  applied  to  first  order  equations  in 
more  accurate  form  by  Heun  [12]  and  Kutta  [14],  and  is 
now  formally  called  the  Runge-Kutta  (RK)  method. 

1.1.1  Basic  Theory  To  solve  numerically  the  initial 
value  problem 

(1.1)  y ’  =  f (x,y) ;  y(x0)  =  yQ 

at  a  series  of  points  x^s  x^  ....  Runge  proposed  the 
following  formula: 


m 


y  ,  ..  -  y  =  T  03 .  k . 
Jn+1  Jn  l  l 


(1.2) 
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where 


(1.3) 


=  hf(xn+aih. 


with  h  =  xn+1  "  xn,  a;L  =  0  and  $  =  0;  j  >  i.  The 

coefficients  on  ,  and  3^  are  chosen  such  that  the 

Taylor  series  expansion  about  (xn*yn)  of  the  on 

the  right  side  of  (1.2)  matches  in  powers  of  h  the 
expansion  of  Yn+1  about  x^;  i.e.  the  left  side  of  (1.2), 


(1.4) 
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m  ^ 

through  terms  in  h  for  some  m*  <  m  (m*=m  for  m<4) 
[21].  In  further  discussion  m*  will  be  referred  to  as 
the  degree  of  the  method,  m  will  be  the  rank  of  the 
method  while  the  term  order  will  refer  to  the  highest 
ordered  derivative  in  the  ODE. 

m  ^  ^  1 

Using  (1.4)  up  to  and  including  terms  in  h  and 

the  differential  operator 


(1.5) 
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where  f  =  f(x  ,y  )  and  j—  ,  j—  are  treated  as  algebraic 

n  n  o  x  o  y 

variables,  an  expansion  for  the  left  side  of  (1.2)  in  terms 
of  f(x,y)  and  its  partial  derivatives  may  be  found.  To 
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obtain  the  expansion  for  the  right  side  of  (1.2)  the 
Taylor  series  for  two  variables,  given  by 


(1.6) 

f ^Xn+hl’yn+h2)  £  i!  B  f 

i  =  0 

where  B  is  the  operator 


(1.7) 

Di  /,  6  ,  ,  6  xi 

B  (hl  Sx  h2  6y  5 

may  be  used.  In  this  case  the  exponent  has  the  meaning 
of  repeated  application  of  B.  To  match  powers  of  h 


through 

m* 

h  it  is  necessary  to  impose  the  restrictions 

(1.8) 

i-1 

ot.  —  £  $,.}  i-2  ,  .  .  .  ,m . 

j=l  J 

1.1.2  RK  Methods  of  Rank  Four  An  example  of  the 
above  procedure  may  be  found  in  Ralston  [ 1 8 ]  for  methods 
of  rank  two,  three  and  four.  After  matching  terms  when 
m  =  m*  =  4  one  obtains  the  constraints  given  in  matrix 
form  by  (1.9). 
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These  constraints  and  the  three  restrictions  imposed 
by  (1.8)  form  a  system  of  eleven  non-linear  equations  in 
thirteen  unknowns.  This  two  parameter  system  can  be 
solved  in  terms  of  and  with  the  obvious  require¬ 

ment  that  the  denominator  not  be  zero  for  any  of  the 
expressions  giving  the  remaining  coefficients.  Several 
methods  of  degree  four  have  been  derived  as  above  simply 
by  choosing  values  of  a 2  and  [19]. 

1.1. 2.1  Truncation  Error  By  maintaining  and 
c 

matching  terms  in  h  Lotkin  [15]  arrived  at  a  bound, 
independent  of  f(x,y) ,  on  the  most  significant  part  of 
the  truncation  error  involved  in  an  RK  method  of  degree 
four.  This  was  accomplished  by  assuming  the  following 
in  a  region  R  about  containing  all  points  in 

(1.3): 


fii+jf,  T,i+j 

(1.10)  |  f  (x,y )  I  <  M;  | -A =H  <~iTT 

6x  oyJ  MJ 

where  i  +  j  <4,  and  L  and  M  are  positive  constants 
independent  of  x  and  y.  The  truncation  error  may  then 
be  written  as 


T, 


-  h5t 


0(h6 ) 


(1.11) 
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where  | t 

(1.12) 


and 

(1.13) 

(1.14) 

(1.15) 

(1.16) 

(1.17) 


is  given  by 

hill  <  [  16  I  u1 1  +  4  I  u2  I  +  |u2+3u3| 

+  |2u2+3u^|  +  |u2+u^|  +  |u^| 

+  8lui|l  +  |u5|  +  |2u5+u?| 

+  |u5+u6+u7|  +  |u6| 

+  | 2u^+u7 |  +  | |  +  2 | Ug | ]ML^ 

ii  h  ii 

=  1/120  -  (a2a)2  +  ag0jg+a^oa^ )/24  , 

2  2 

u2  _  1/20  ~  [ a2a g$ g2°°3  i_j. ^ 2^42 3^43^ 4 

3  3  3 

Ug  =  1/120  -  [a23g2c0g  +  (a23ij2+ag6i4g)a)^]/6 

2  2  2 
=  1/30  -  [a2oig3g2a)g  +  (a23i|2+a3^43^c04 

Ug  =  1/120  -  a23g231|g001]  , 

u6  =  1/40  -  [a23g2a)3  +  (a231|2  +  a3343)2w4]/2 


]/2  , 


]/2  , 


(1.18) 
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(1-19)  =  7/120  -  a2(a3+ai4)6326i|3wi|  , 


(1-20)  Ug  =  1/120  . 


Ralston  [ 1 8 ]  proceeded  further,  obtaining  coefficients 
for  methods  of  degree  two,  three  and  four  which  minimized 
this  type  of  bound. 

1.1.3  Advantages  and  Disadvantages  Ralston  [19] 
cites  the  advantages  and  disadvantages  of  RK  methods. 

The  three  most  advantageous  properties  are : 

1.  They  are  single  step  methods,  i.e.  to  find 
yn+1  it  is  necessary  only  to  have  the  pre- 
ceeding  values  of  yn  and  x  and  the  step- 
size  h. 

2.  They  are  self-starting,  the  interval  between 
steps  may  be  changed  at  will  and  they  are 
particularly  easy  to  apply  on  a  digital  computer. 

3.  They  are  comparable  in  accuracy  to  predictor 
corrector  methods  of  the  same  degree. 

The  drawbacks  to  these  methods  are  twofold: 

1.  The  step-size  generally  must  be  chosen  conser¬ 
vatively  because  of  the  difficulty  in  obtaining 
the  per-step  error,  i.e.  h  is  smaller  than  is 
actually  necessary  in  order  to  achieve  the 
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desired  accuracy. 

2.  They  require  a  number  of  evaluations  of  f(x,y) 
at  each  step  equal  to  or  greater  than  the  degree 
of  the  method,  whereas  the  corresponding  pre¬ 
dictor-corrector  method  requires  a  number  of 
function  evaluations  less  than  the  degree. 

It  is  because  of  the  above  restrictions  that  RK 
methods  are  in  general  used  only  to  obtain  starting 
values  for  predictor-corrector  methods  and  in  changing 
the  step-size. 

1 . 2  Summary  of  Literature 

Work  has  continued  on  RK  methods  for  first,  second 
and  higher  order  equations.  Bieberbach  [1],  as  well  as 
Lotkin,  did  extensive  research  on  the  truncation  error 
involved  in  these  methods  while  Blum  [2],  Ceschino  and 
Kuntzmann  [3],  and  Gill  [10]  presented  new  sets  of  coef¬ 
ficients  having  specific  attributes.  Nystrom  [16] 
developed  methods  of  degree  greater  than  four  and  Frey 
[7],  who  refers  to  these  higher  degree  methods  as  Runge- 
Kutta-Nystrom  methods ,  improved  upon  them  by  choosing 
coefficients  such  that  the  number  of  function  evaluations 
necessary  for  each  step  is  less  than  the  rank  of  the 


method  for  m  >  5 . 
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Nystrom  [ 1 6 ]  also  extended  RK  methods  to  second 
order  equations  and  Zurmuhl  continued  this  extension, 
first  to  third  order  equations  [22],  and  then  to  equa¬ 
tions  of  the  n-th  order  [23].  Nystrom  restricted  his 
work  for  second  order  equations  to  methods  of  rank  three 
and  four  while  Zurmuhl  concerned  himself  only  with  methods 

of  rank  four  and  function  evaluations  at  x  ,  x  +  h/2 

n s  n 

and  x^  +  h.  Nystrom  [16]  provides  a  set  of  coefficients 
for  second  order  equations  while  Zurmuhl  [23]  presents 
computing  schemes  for  n-th  order  equations  where 
3  *  n  <  10. 

Froese  [9]  evaluated  the  RK  type  methods  suggested 
by  Zurmuhl  using  Rademacher’s  theory  of  error  propaga¬ 
tion  [17].  Her  work  was  primarily  concerned  with  the 
fourth  order  initial  value  problem 

(1.21)  y^  =  y;  y^(xQ)  =  y^;  i=0,...,3, 

and  the  question  of  whether  it  should  be  solved  as  a 
system  of  four  first  order  equations,  two  second  order 
equations  or  by  use  of  a  direct  method  for  fourth  order 
equations.  Zurmuhl  [23]  states  that  solving  an  n-th 
order  equation  directly  is  superior  to  solving  an  equiva¬ 
lent  system  of  equations  if  n  >  3.  This  conclusion  was 
based  on  a  consideration  that  the  truncation  error  for 
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the  RK  method  of  rank  four  is  of  degree  (n+2)  for 
n  >  2.  Froese’s  analysis  and  experimental  results  con¬ 
firmed  Zurmuhl’s  conclusion. 

1. 3  Purpose  of  Study 

Any  direct  method  for  second  order  equations  will 
hereafter  be  called  a  Runge-Kut ta-Ny strom  (RKN)  method. 
The  purpose  of  the  present  investigation  is  to  present 
a  new  derivation  of  the  general  RKN  method  of  rank  m 
and  degree  m* »  and  to  obtain  several  new  sets  of  coef¬ 
ficients  for  methods  of  rank  four.  The  coefficient 
constraints  for  RKN  methods  of  rank  three  and  four  are 
derived  with  a  new  restriction  reducing  the  number  of 
equations  in  these  systems.  Bounds  for  the  truncation 
error  of  these  methods  are  also  found  using  the  approach 
of  Lotkin  [15]. 

Several  second  order  ODE’s  are  solved  using  a  system 
of  two  first  order  equations  and  the  corresponding  RKN 
methods.  Comparing  these  results  and  the  bounds  on  the 
truncation  error  involved  in  both  methods  it  is  shown 
that  Zurmuhl's  result  [23]  could  be  extended  to  include 
n  =  2 . 
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CHAPTER  II 


RUNGE-KUTTA  METHODS  FOR  SOLUTION 
OF  SECOND  ORDER  ODE’s 

2 . 1  As  a  System  of  First  Order  Equations 
An  r-th  order  ODE 

(2.1)  =  f  (x,y  ,y  ’  ,  .  .  .  ,y  (r_1^  ) 

y(l)(x0)  =  i  1  =  0,1,  ...  ,r-l 

may  be  solved  as  a  system  of  lower  order  equations  if  so 
desired  [9].  For  simplicity  the  order  of  all  the  equa¬ 
tions  in  the  system  may  be  chosen  to  be  one.  In  this 
case  (2.1)  becomes 


(2.2) 


Y-  =  F(x , Y) ;  Y(x0)  -  YQ 


and  may  be  solved  by  applying  RK  methods  [11] 


(2.3) 


m 

Y  , -,  =  Y  +  I  u).K, 
n+1  n  l  l 


where 


K.  =  hF(xn+a.h, 


(2.4) 


o.l 


! 


’  . 


i  . 


12 


The  variables  Y 

r ,  i  .  e  .  K.  = 

1 

Hildebrand 
value  problem 


K 


and  K.  denote  vectors  of  length 


( kH  ?ki2  » •  •  •  >kir)  • 

[13]  considers  the  second  order  initial 


(2.5)  y"  =  f(x,y,yf);  y(xQ)  =  yQ;  y’(xQ)  =  y^  , 
treating  it  as  the  system  of  first  order  equations 

(2-6)  yT  =  z;  y(x0)  =  y0 


and 


(2.7) 


z'  =  f(x,y,z);  z(xQ)  =  zQ  , 


and  suggesting  the.  solution  by  means  of  a  RK  method  of 
rank  four.  Extending  this  to  a  method  of  rank  m  would 
give  the  solution  of  (2.6-7)  and  consequently  (2.5)  in 
the  following  manner: 


(2.8) 


y 


n+l 


m 

y  +  T  u.k. 
i=l 


and 


z 


n+l 


m 

z  +  J  a).  £ . 
i=l 


(2.9) 


■ 


. 


' 


. 
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where 


(2.10) 

ki  = 

and 


(2.11) 

£.  =  hf(xn+a1h,  yn+.i  Ajkj>  V  .i  AjV  • 

J  J 

The  above  computational  procedure  where  m  =  4  is  that 
most  commonly  used  for  solving  (2.5)  using  RK  methods. 

2 . 2  Derivation  of  the  RKN  Method 

The  remainder  of  this  Chapter  presents  a  new  deriva¬ 
tion  of  the  general  RKN  method.  Nystrom  [ 1 6 ]  derived  RKN 
methods  of  rank  three  and  four  while  Zurmuhl  [22],  [23] 
worked  only  with  methods  of  rank  four  for  third  and  higher 


ordered 

equations.  To  develop  the  general  RKN  method  of 

rank  m 

define 

(2.12) 

v  0 , 

v  ’  =  —  •  i  =  —  k* 
yn  h  5  i  hi' 

Substituting  (2.12)  into  (2.10)  gives 


(2.13) 

i-1 

k.  =  v  +  2!  7  6. .k*  , 
i  n  A  ij  J  ’ 

(2.13) 


9 


- 


. 


-  :'i 


' 


* 


o 
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and  subsequently  using  (2.12-13)  in  (2.11)  yields 


(2.14)  i.  =  21 


i-1 


i  h  k!  ■  hf(*n+oih*  yn+( 


j=i  -  p 

From  (2.13-14)  (2.8-9)  becomes 


j=l 


m  mi-1 

(2,15)  yn+l  =  yn  +  (  l  "i)vn  +  2!  1  "l (  I  6liki) 

i=l  n  i=l  1  j  =  l  J 


and 


(2.16) 


z  , =  z  +  -T-  J  oi.k* 
n+1  n  h  . L n  i  i 

i=l 


m 


Imposing  the  following  constraints  included  in  any  RK 
method  of  rank  m: 


i-1 

(2. 

17) 

ou  =  I 

J=1 

and 

(2. 

OO 

i — 1 

i  =  2  5  3, . . • ,m 


m 


y  a).  =  i  . 

l  i  ^ 


1=1 


and  multiplying  (2.16)  by  h  gives  for  (2.15-16) 


' 


'  ;  '  '■ 


N  • 


• 

. 
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(2.20) 


V 


n+1 


v 

n 


m 


+  2!  V  co.  k* 
L  11 
1=1  x  1 


where 


h2  1-1  J-1 

(2.21)  k*  =  -^-r  f(x  +a.h,  y  +a.v  +2!  7  3..(  I  3.  k*), 
l  2!  n  i  5  Jn  l  n  ^  ij  ^  OP  P  5 


i-1 

(v^+2!  £  3, ,k?)/h)  . 

j 


n  ij  J 


Applying  (2.19-21)  directly  to  (2.5)  Is  thus  equiva¬ 
lent  to  using  (2.8-11)  in  solving  (2.6-7).  The  above 
leads  to  the  formulation  of  the  general  RKN  method  which 
may  be  conveniently  given  by 

m 

(2.22)  y  , -,  =  y  +  v  +  I  wn.k. 

v  J  n+1  J  n  n  1,  Oi  l 

i=l 


and 


(2.23) 


v 


n+1 


m 


T/ii 

i  =  l 


where 


■ 


1  . 
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(2.24) 


1-1 


(v  +2!  7  c. .k. )/h) . 
n  .4,  ij  j 


2.2.1  Coefficients  for  an  RKN  Method  Matching 
corresponding  terms  in  (2.19-21)  and  the  general  RKN 
method  given  by  (2.22-24)  yields  coefficients  for  a 
RKN  method  of  rank  m  in  terms  of  the  coefficients  of 
any  RK  method  of  corresponding  rank.  These  are  given  by 


0;  i=l,2,...,m;  j>i-l 


i-1 


(2.25) 


cij  1  j  1>2 


m 
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These  relationships  between  RK  and  RKN  methods  enable 
one  to  obtain  new  coefficients  for  RKN  methods  of  any 
rank  if  the  coefficients  for  an  RK  method  of  the  same 
rank  are  known. 

2.2.2  Coefficient  Format  The  coefficients  for 
RK  methods  of  rank  m  will  be  given  throughout  in  the 
following  form: 


(2.26) 


The  corresponding  RKN  format  will  be 


3.i?  f  ir J 


- 1 


18 


(2.27) 


b21 

2 

a^ 

i — i 

on 

rQ 

b  32 

2 

a3 

0 

• 

• 

bmi 

• 

• 

• 

bm2 

•  •  • 

b  , 

m,m-l 

• 

• 

• 

2 

a 

m 

■ — i 

OJ 

o 

a  ^ 

1 — 1 
on 
o 

c32 

• 

a3 

• 

i — i 

•  •  £ 

o 

• 

Cm2 

•  •  • 

m  ,m-l 

• 

a 

m 

W01 

W0  2 

o  •  • 

Wn  i 

0  ,m-l 

Wn  rr, 

u  ,m 

i — 1 
i — 1 
& 

W12 

•  e  • 

W-,  i 

1  ,m-l 

W-, 

1 3m 

The  reason  for  using  this  format  will  become  clear  in 
the  following  Chapter.  All  coefficients  restricted  to 
a  value  of  zero  are  not  included. 


■ 


CHAPTER  III 


DERIVATION  OF  CONSTRAINTS  FOR 
RUNGE-KUTTA-NYSTROM  METHODS 

3 • 1  Constraints  for  a  Method  of  Rank  Four 

The  classic  derivation  of  the  RK  method  of  rank 
four  was  outlined  briefly  in  Chapter  I  to  indicate  the 
manner  in  which  the  technique  was  originally  formulated. 
The  same  pattern  will  be  followed  in  more  detail  in 
building  similar  constraints  for  RKN  methods  of  the  same 
rank . 

Consider  (2.22-23)  expressed  in  the  form 

m 

(3.1)  y  ,  -  y  -  v  =  J  wn.k. 

v  J  '  Jn+1  Jn  n  Oi  i 

and 

m 

(3.2)  v  -  v  =  I  w,.k.  . 

v  '  n+1  n  li  i 

To  derive  the  constraints  for  a  RKN  method  of  rank  m 
it  is  necessary  to  expand  the  left  and  right  sides  of 
(3.1-2)  by  means  of  Taylor  series.  The  coefficients  of 
the  method  are  chosen  in  such  a  manner  that  the  expan¬ 
sions  of  each  side  of  (3-D  Have  matching  terms  in  powers 

Q  rn  HI  ^  1 

of  h  through  h  and  in  (3.2)  through  h  ,  the 


, 

_ 

H  \  d 


■ 
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latter  because  of  the  definition  for  given  by  (2.12). 

3.1.1  The  Taylor  Series  Expansion  Expanding  yn+p 
about  x  =  using  the  Taylor  series  for  a  function  of 
one  variable  given  by  (1.4)  and  replacing  hy^  with  vn 
according  to  (2.12)  produces 


(3.3) 


y 


n+1 


-  y 


n 


v  = 
n 


h_ 

2! 


y 


n 


h_ 

3! 


y 


(3) 


n 


h4 

TT 


(4) 


n 


h' 


5!  yn 


(5) 


for  the  expansion  of  the  left  side  of  (3.1).  Similarly 
the  Taylor  series  expansion  of  vn+p  about  x  =  xn 
with  derivatives  of  vn  replaced  by  the  derivatives  of 
yn  from  (2.12)  gives 


(3.4) 


v 


n+1 


-  v 


n 


=  h2y" 
j  n 


h_ 

2  ! 


3 


y 


(3) 


n 


+ 


h_ 

3 


y 


(4) 


n 


+ 


(5) 

n 


+ 


0  0  0 


for  the  left  side  of  (3.2). 

( 3 ) 

For  the  method  of  rank  four  expressions  for  y 
y(^  and  y^  must  be  found  in  terms  of  f(x,y,y'). 
For  compactness  one  may  use  the  differential  operator 


' 


'  I 
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(3.5) 

-n1  /  6  ,  6  ,  n  (5  J 

D  =  (t —  +  z-7 —  +  f*? — ) 

6x  6y  6z 

where 

6  6  6 

z  =  y'  f  =  f (x  ,y  .y’)  and  -r— ,  and  -j-r 

°n  n,Jn,Jn  6x5  6y  <5z 

are  treated  as  algebraic  variables.  Thus 


(3.6) 

=  A  f  =  f  +  Zf  +  ff  =  Df  , 

J  n  dx  x  y  z  3 

(3.7) 

y^'J  =  Df  =  f  +  2zf  +  2 f  f 

Jn  dx  xx  xy  xz 

+  2zff  +  z2f  +  f2f 

yz  yy  zz 

+  f  (f  +zf  +ff  )  +  ffT 
z  x  y  z  y 

=  D2f  +  f  Df  +  ff 

z  y 

and 


(3.8) 

y(5)  =  d  (d2  +  f  Df  +  ff.  ) 

dx  z  y 

f xxx  +  3zfxxy  +  3ffxxz 

+  z  fyyy  +  3z  fxyy  +  3z  ffyyz 

+  6zffxyz  +  3f2fxzz  +  3zf2fy^ 

( „  „  .  +  ■ 


. 


' 


+  f°f  +  3f(f  +zf  +ff  ) 
zzz  xy  yy  yz' 


+  3(f  +zf  +ff  )(f  +zf  +ff  ) 
x  y  z  xz  yz  zz 


+  f  „[f  +  2zf  +  2ff 

zL  xx  xy  xz 


+  2zff  +  z2f  +  f2f 

yz  yy  zz 


+  f  (f  +zf . +ff  )  +  ff  ] 
z  x  y  z'  y 


+  f  (f  +zf . +ff  ) 
y  x  y  z' 


-  n3 


D  f  +  3fDf  +  3DfDf 
y  z 


where 


+  f  (D  f+f  Df+ff  )  +  f  Df 
z  z  y  y 


(3.9) 


6nf 


tn  t0  .  .  .  t  <5t-,  fit 0  .  .  .  fit 
12  n  12  n 


The  substitution  of  (3.6-8)  into  (3*3-*0 


(3.10) 


2  3 

h  hJ 

y  -  y  -  V  =  TT-r  f  +  rr-r-  Df 

Jn+1  Jn  n  2!  3! 


,4  p 

+  (D^f+fzDf+ff  ) 


gives 


’  • 


-  ■ 


23 


+  j=y  [D3f  +  3fDf  +  3DfDf 

^  y  ^ 

+  f  (D^f+f  Df+ff  )  +  f  Df] 

z  z  y  y 

+  0(h6 ) 


and 


3  M 

(3'11)  vn+l  -  vn  =  h2  f  +  b  Df  +  JT  (D2f+fzDf+ffy) 

h5  , 

+  T^j-  [D  f  +  3fDf  +  3DfDfz 

+  f  (D2f+f  Df+ff  )  +  f  Df] 
z  z  y  y 

+  0(h6) 

which  are  the  required  expansions  for  the  left  sides  of 
(3.1-2)  in  terms  of  D. 

3.1.2  Expansion  of  the  To  obtain  the  expan¬ 

sions  for  the  right  sides  of  (3.1-2)  the  Taylor  series 
for  a  function  of  three  variables  about  the  point 

(x  .y  .y’)  is  used.  It  is  given  by 
n3Jn>Jn  ° 

m#  i  ^ 

(3.12)  f(xn+hl»yn+h2’yA+h3)  =  I0  IT  B  f  +  0(h  5 


- 

L 


. 


' 


. 
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where 


(3.13) 

•pi  _  /u  6  iv,  5  .  ,  6\i 

(hl  6x  h2  6y  +  h3  Iz)  > 

z  =  as  defined  previously  and  the  exponent  has  the 

meaning  of  repeated  application  of  B.  The  k±  * s  for 
a  method  of  rank  four  are 


(3.1*0 

kl  "  2!  f(xn,yn,yn)  =  5T  f  » 

(3.15) 

h2 

k2  =  2T  f[xnta2h-  yn+a2vn+b21kl>  { V2o21k1)/h]  > 

(3.16) 

h2 

k3  =  21  f[Va3h>  yn+a3vn+b31kl+b32k2* 

(vn+2(c31k1+c32k2) )/h] 

and 


(3.17) 

h2 

k 4  -  2  !  yn+a4vn+b4lkl+b42k2+b43k3 > 

( vn+2(cZllk1+ci|2k2+cij3k3)  )/h]  . 

Since  a  method  of  rank  four  is  also  of  degree  four 

2 

it  is  necessary,  in  view  of  the  factor  h  /2!,  to  retain 

1  2 

only  terms  through  h  when  expanding  2!k^/h  .  Sub- 


' 


.  :  '  .  i  ;  ric  b  3J£-  >r3  3  •  ax. s  ■ 

■ 


. 


stituting  (2.12)  and  (3.1*0  into  (3.15)  and  expanding 
yields 
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(3.18) 


2  *  Vi 

k2  =  f[xn+a2h,  yn+a2hz+fjb21f »  z+c2ihf^ 


=  f  +  [a2hfx  +  a2hzfy  + 


+  c21hffZ^  +  FT  C<a2h)2fxx 


+  (a2hz)2fyy  +  2(a2hz)(|Tb21f)f 


yy 


+  (c2]hf)  fzz  +  2(a2h) (a2hz)fxy 


-h 


+  2(a2h)(^yb21f)fxy  +  2(a2h)(c21hf)f 


xz 


h‘ 


+  2(a2hz) (c21hf )fyz  +  2 (^jb21f ) (c21hf ) fyz ] 


yz 


+  IT  [(a2h)3fxxx  +  (a2hz)3fyyy  +  ( °21hf > 3f zzz 


+  3(a2h)2(a2hz)fxxy  +  3(a2h) (a2hz)2fxyy 


+  3(a2h)2(c21hf)fxxz  +  3(a2h)(o21hf)2fxzz 


+  3(a2hz)2(o21hf)fyyz  +  3(a2hz)(c21hf)2fyzz 


+  6(a2h) (a2hz) (c21hf)fxyz]  . 


' 


'  ■ 


: 


:  •  o:  --  - 


:  ■ 


■ 


. 
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h 


Gathering  terms  in  jy  ;  i=0,l,2,3  gives 


2  ! 


(3.19)  -j  k2  =  f  +  h[a2fx  +  a2zfy  +  c2iffz-* 


+  IT  Ca2fxx  +  4*\y  +  c21f2f 


zz 


+  2a2zfxy  +  2a2c21ffxz  +  2a2c21zffyz 


+  b21ffy]  +  IT  [alfxxx  +  aIz3fyyy 


+  °21f3fzzz  +  3aizfxxy  +  3aiz2fxyy 


+  3a2c21ffxxz  +  3a2o21ffxzz 


+  3a2c21z2ffyyz  +  3a2c21zf2fyzz 


+  6a2o21zffxyz  +  3a2b21ffxy 


+  3a2b21zffyy  +  3b21c21f  fyz^  * 


Utilizing  the  operator  D  as  given  by  (3.5)  and  the 
restrictions 


a2  =  c21;  b21  =  a2 


(3.20) 


W  ■  f  m”  b  If  +  Ynisd  + 


is  VE  it.  - ;°<cx 


QClOlJOl  3  'I 
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one  obtains 


(3.21) 


k2  =  b  f  +  It  a2Df  +  2TTT  a|(D2f+ff  ) 


+  573T  a3(D3f+3fDf  )  +  0(h6) 


Substitution  of  (3. 11*)  and  (3.21)  into  (3.16)  re¬ 


taining  terms  through  h5  gives 


(3.22) 


2  |  ^ 

2  —  *  yn+a3bz+2Tb31f 


2  3 

+  b32(§Tf+§Ta2Df) *  z+c31hf 


+  c32(hf+a2h2Df+|ya2(D2f+ffy))] 


h‘ 


=  f[xn+a3h,  yn+a3hz+(b31+b32)^yf 


+  ^-j-a2b32Df,  z+(c31+c32)hf 


+a3c32h2Df+|ya2c32(D2f+ffy))]  . 


Expanding  in  the  same  manner  as  2!k2/li  and  imposing 
the  restrictions 


(3.23) 


a3  —  c32 »  —  b3i  b32 


'  •  '  Q  - 


. 

■ 

■ 


. 


•xeariST  i  ei  .  &•*.  J  1  j  c. 
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gives  in 

terms  of  D 

(3.24) 

2  3  4 

k3  =  Jr  f  +  |t  a3Df  +  2J2T  Ha|(D2f+ffy) 

+  2a2c32fzDf]  +  gT^T  [a3(D3f+3fDf  ) 

+  3a2C32(fzD2f+ffyfz)  +  6a2a3c32DfDfz 

+  3a2b32fyDf]  . 

Following  a  procedure  similar  to  that  for  k2  and 


k^  with 

the  restrictions 

(3.25) 

2 

a4  ”  c4i  +  c42  +  c43’  a2  =  b4l  +  b42  +  b43 

yields 


(3.26) 

k4  =  2!  f  +  2!  a4Df  +  5T2T  C^(D2f+ffy) 

5 

+  2 ( a^c^2+a3c^3 ) f zDf ]  +  2! 3!  [&^(D  f+3fbfy) 

+  3(a2bi)2+a3b1)3)fyDf  +  3(a2o1)2+a2o1J3)fz(D2f+ffy  ) 

2 

+  6a^ ( a2c ||2+a^c ^ ) Df Df z  +  bd.2^^2^  • 

' 


1 


, 


£ 


■ 


■ 


■  :C  .  £  +  -  (l  :  •  ;  ;r. 
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Multiplying  (3.14),  (3.21),  (3.24)  and  (3.26), 
respectively,  by  the  corresponding  wpl,  p=0,l,  and 
gathering  terms  gives  the  following,  respective, 
expansions  for  the  right  sides  of  (3.1-2): 


(3-27)  wplkl  +  Vk2  +  wp3k3  +  wp4kH 


=  (wpl+wp2+wp3+wpi))  2 T  f 


+  (a2Wp2+a3Wp3+a1)Wp1))  jy  ™ 

4 

+  (a2wp2+a3wp3+a4wp4^  2l2T  f+ffy^ 

h4 

+  2[a2c32wp3  +  (a2ci|2+a3cij3)wp4]  fzDf 

5 

+  ^a2wp2+a3wp3+a4wp4^  7TTT  (C)3f+3fDf'y) 

f 

h^ 

+  3[a2c32wp3  +  (a2ci|2+a3clj3)wpi|]  JTJT  fz^D 


h' 


+  3Ca2b32wp3  +  (a2b42+a3b43^Wp4^ 2 ! 3!  fyDf 


+  6[a3a2c32wp3  +  a„ (a2ci|2+a3c1(3)wpi)]  2737 


h5  2 

+  6a2o32o1)3wplt  ^fjT  fzDf 


2f+ff  ) 


DfDfz 


. 

!t3  r  (|lqw(lB+£qwEB+SqwSB> 

• 

' 

• 
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where  p  =  0,  and 


(3.28)  same  as  (3.27)  except  p  =  1 

2  5 

Matching  terms  in  powers  of  h  through  h 
(3.10)  and  (3.27),  and  (3.11)  and  (3.28)  yields  the 
constraints  given  in  matrix  form  by  (3.31).  The 
restrictions  used  in  finding  the  expansions  for  the 
are  necessarily  included  as  constraints  and  are 


given  by 

i-1 

(3.29) 

ai  =  ^  c±]; 

j=l  J 

i‘2,3,4, 

and 

2  i-1 

(3.30) 

'  ai  =  £  bi1» 

1  J-l  J 

1=2, 3, 

for 


wS  ab  1  iXx  (  !S  .£)  ons  (I 
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1 - 

OJ 

on 

-=r 

lo 

LTv 

LT\ 

LT\ 

LT\ 

— \ 

1 — 1 

1 — 1 

1 — 1 

1 — 1 

» 

i — 1 

1 - 1 

1 - 1 

rH 

rH 

1 _ 

i 

ii 

1 

i 

i — 1 

CM 

on 

-=r 

i — 1 

1 - 1 

rH 

l — 1 

5 

£ 

5 

rH 

CM 

on 

-=r 

o 

O 

o 

o 

£ 

1 _ 

£ 

5 

j 

1 

1 

on 

-=r 

on 

on 

on 
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Other  derivations  for  constraints  of  the  RKN  method 
of  rank  four,  notably  Nystrom  [16]  and  Duffin  [5],  ex¬ 
cluded  the  restriction  imposed  by  (3.30).  In  this 
instance  six  other  constraints  are  found  as  follows: 

h4 

from  Tpr  ff 
4 !  y 


(3.32) 

6Cb21w02  + 

c\j 

on 

X> 

+ 

rH 

on 

)w03  + 

(b4^+b42+b43 )Wq4 ]  =  1 

(3.33) 

6Cb21w12  + 

(b  3i+b  32 

)w13  + 

(b 4 1+^ 42"*"^ 4 3 ^ wl4 ^  _  ^ 

h5 

from  ft  fDf 

b !  y 

( 3.  3*0  30 [a2b21WQ2  +  a^(fr 31"**^ 32^w0 3  a4^4l"^^>42"^^)43^w04^  -  ^  * 

(3.35)  30[a2b21w12  +  a3(b31+b32)w13  +  ai4(b1)1+b1|2+bi|3)wl4]  =  5  , 

h5 

from  ^  ffy 

(3.36)  30[b21c32wQ3  +  (b21o1)2+(b31+b32)c1)3)]w0l)  =  1  , 

(3.  37)  30[b21o32w13  +  (b21o1)2+(b31+b32)o1)3)]w11)  =  5  , 

and  would  have  been  included  in  (3.31).  As  is  easily 
seen  the  imposition  of  (3.30)  reduces  the  number  of 


v 


" 


■ 
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constraints  involved  in  a  method  of  rank  four  from  twenty- 
seven  to  twenty-four. 


3 . 2  Truncation  Error  for  a  Method  of  Rank  Four 

The  truncation  error  involved  in  a  RKN  method  of 
rank  m  is  given  in  general  by 


(3.38) 


hm*+1 

^m  ”  (m*+l ) !  ^m 


0(hm*+2) 


The  procedure  followed  by  Lotkin  [151  and  outlined  in 
Chapter  I  for  a  RK  method  of  rank  four  will  be  used  in 
obtaining  a  similar  bound  on  the  truncation  error  for  a 
RKN  method  of  the  same  rank. 

All  constraints  included  in  (3.29-31)  with  the 
exception  of  those  given  by 


(3.39) 


0 

10a2 

ioa! 

lOajj 

l 

1 — 1 

O 

£ 

1 

1 

0 

0 

2 

3  0  3*  ^  ^  32 

2  2 

30  ( a  2^  2"t  a  ^ 

W02 

1 

0 

0 

30a2b^2 

30  (a2b/|2+a3b^2 ) 

W03 

1 

0 

0 

60ct^a.2^'^2 

60a^  (a^c^^ta^c^i^) 

w0i4 

1 

0 

0 

0 

60a2^32^i|3 

1 

. 


. 
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must  be  satisfied  for  a  RKN  method  of  rank  four.  It  is 

these  constraints  which  determine  the  value  of  t  in 

m 

(3.38)  for  m  =  m#  =  4 . 

Throughout  the  derivation  of  the  constraints  in 

5 

Section  3.1  terms  in  h  were  retained,  not  only  because 
necessary  constraints  were  derived  from  them,  but  also  to 
obtain  a  bound  on  |tjJ  .  To  obtain  an  estimate  of  the 
truncation  error  involved  it  is  necessary  to  find  the 
difference  between  terms  in  h^/5!  in  (3.10)  and  (3-27). 
This  is  given  by 


(3.40)  t4  =  u1(D3f+3fDf  )  +  u2fz(D2f+ffy) 

+  u3fyDf  +  u4DfDfz  +  u5f2Df  , 

where 


o. 41) 

(3.42) 

(3.43) 

(3.44) 


3  3  3 

u±  =  1  -  10(a2W02+a3w03+ai|w0^)  , 

2  2  2 

u 2  =  1  —  30(^2^22^02  ( ^2 ^ 4 2"^^3^ 4 3 ^0 4 ^  * 

u^  =  1  -  30 [a2b 22w03  +  (a2b42+a3b43^w04-^  5 
u^  =  1  -  60[a3a2c32w03  +  a^ ^c^+a^^w^] 


and 


•y  s e.  9  '  .  u*i : 
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(3.45)  Up-  —  1  —  60  (a2c^2c^^w0^ )  . 

The  evaluation  of  by  (3.40)  is  dependent  upon 

?(x>y>y')  while  the  purpose  of  Lotkin's  procedure  is  to 
obtain  a  bound  on  the  truncation  error  independent  of 
this  function.  This  is  accomplished  by  defining  the 
following  in  a  region  H  containing  all  points  in  (2.24) 


(3.46) 


I f (x,y ,y ’ ) |  <  M;  | y 1 |  =  | z |  <  N 


^i+j+kf 
5x^6y^  6zk 


Li+j+k 

nV-1 


where  i  +  j  +  k  <  4  and  L,  M  and  N  are  positive 
constants  independent  of  x  and  y. 

Using  (3.46)  enables  limits  to  be  placed  on  the 
terms  in  (3.40)  dependent  on  f(x,y,yf).  The  basic 
elements  in  these  terms  are  given  below  with  their 
respective  limits. 


(3.47) 
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f 

f 
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f  <  l2m 
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2  -1 
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Cm 
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lV1 
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XZ 

yz 

<  l3mn_1 

f  <  l3mn"2 

f  vvr7 

no 
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xxy 

xyy 

xxz 

3  -1 
<  L3M 

3  -2 

f  <  LrN 

f 

<  L3M  1N 

xzz 

yyz 

y  zz 

r.o  \  n  r- 
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The  seven  coefficients  of  the  on  the  right  side  of 

(3.^0)  may  be  rewritten  without  D  as 


(3. ^8) 

D3f  =  f  +  z3f  +  f3f  +  3zf 

xxx  yyy  zzz  xxy 

+  3ff  +  3z2f  +  3f2f  +  3z2ff 

xxz  xyy  xzz  yyz 

+  3zf2f  +  6zf f  , 

yzz  xyz  * 

(3.^9) 

fDf  =  ff  +  z  f f  +  f2f  , 

y  xy  yy  yz  3 

(3.50) 

fzc2f  ■  Vxx  +  zVyy  +  f2fzfzz 

+  2zfzfxy  +  2ffzfxz  +  2zffzfyz  ' 

(3.51) 

ff  f  =,ff  f  , 
y  z  y  z 

(3.52) 

f  Df  =  f  f  +  zf2  +  ff  f  . 
y  x  y  y  y  z  3 

(3.53)  DfDfz  =  fxfxz  +  zfxfyz  +  ffxfzz  +  zfyfxz 


■  '  1  'V  : 
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and 


(3.54) 


fzDf  “  fzfx  +  zfzfy  +  ffz 


Collecting  terms  on  the  right  sides  of  (3.48-5*0  and 
using  their  respective  limits,  easily  obtained  from 
(3.47),  the  following  absolute  bound  on  (3.40)  is  ob¬ 
tained  : 

(3.55)  | t4|  <  (27|u1|  +  4 | u2 |  +  |u2+u4|  +  2|2u2+u4|  +  6  |  U||  | 


+  3|u5|)L3M  +  (9|u1|  +  |u2+u3|  +  2 | u^ | ) L^M^N 


2,.2xt-1 


For  a  RKN  method  of  rank  four  the  absolute  value  of  (3.38) 
can  then  be  written  as 


2..2xt-1 


(3.56)  |T^|  <  h"5  (e1L'5M  +  e2L  M  N”  )  +  0(h  ) 


where 


(3.57)  e±  =  (27|u1|+4|u2|  +  |u2+u1||+2|2u2+ui|| 


+6|u^|+3|u5| )/5 ! 


and 


. 
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(3.58)  -  ( 9  |  u-^  |  +  |  U2+u^  | +2  |  |  ) /5  !  . 

A  second  order  ODE  often  encountered  is  that  given  by 

(3.59)  y”  =  f(x,y);  y'(xQ)  =  y^;  y(xQ)  =  yo 

in  which  the  first  derivative  y*  does  not  appear 
explicitly.  In  this  case,  since  any  partial  derivative 
with  respect  to  z  is  equal  to  zero,  the  expression  for 
the  truncation  error  involved  in  solving  (3.59)  by  a  RKN 
method  is  significantly  different  than  that  given  by 
(3.56-58).  If  all  terms  are  ignored  which  contain  a 
partial  derivative  with  respect  to  z  it  is  found  that 

(3.60)  1 1 4 |  <  8|u1|L3M  +  (6|u1|+2|u3| )L2M2N_1 
and 

b5  o 

(3.61)  |TJ  <  [8 1  u  |L% 

f (x,y )  i 

+  (6|u1|+2|u3| )]L2M2N_1  +  0(h6 )  . 

3. 3  The  RKN  Method  of  Rank  Three 

Using  the  expansions  for  the  left  and  right  sides 
of  (3.1-2)  derived  in  Section  3.1  for  m  =  4  the 
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corresponding  expansions  for  the  RKN  method  of  rank 
three  can  be  obtained  simply  by  retaining  only  terms 

li 

through  h  .  Prom  (3.10-11) 


(3.62) 


h2  h3 

yn+l  -  yn  '  vn  =  2T f  +  JTDf 


4 

+  ^j(D2f+fzDf+ff  )  +  0(h5) 


and 


3  4 

(3.63)  vn+1  -  vn  -  h2f  +  lyDf  +  |T(D2f+fzDf+ffy)  + 


while  (3.27-28)  gives 


(3.64-65) 


w  ,k-.  +  w  0k0  +  w^0k_ 
pi  1  p2  2  p3  3 


"  (wpl+wp2+wp3)  lrf 


h 


3 


+  (a2wp2+a3wp3^  2TDf 


+  (a2wp2+a2wp3)  2l2T(D2f+ffy) 


+  2a~cnow 


2  32  p3  2 ! 2 !  z 


f  _Df 


0(h5)  , 


where  again  p  =  0,1  . 
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3.3.1  The  Constraints  for  a  RKN  Method  of  Rank 

2  4 

Three  Matching  terms  in  powers  of  h  through  h 
for  (3.62)  and  (3.64),  and  (3.63)  and  (3.65)  yields  the 
necessary  constraints  given  in  matrix  form  by 


& 
2  ! 1 


hiDf 

3! 


(3.66) 


h4  5 

^jy(D  r+ff  ) 


TTfzDf 


1  1 


0  3a2  3a3 


2  2 
0  6a2  6a^ 


0  0 


12a^ ^^2 


W01  W11 


W02  W12 


wQ3  w13 


1  2 


1  3 


The  necessary  restrictions 


(3.67) 


i-1 

=  I  ci-<  >  1  =  2 ,3 

J-l  J 


and 

9  i-i 

(3.68)  af  =  l  b  ;  i=2,3 

1  j-l  10 

must  be  included  with  (3.66)  to  comprise  the  full  set  of 
constraints  for  a  method  of  rank  three.  If  (3.68)  is 
not  imposed  two  other  constraints  are  found  from  the  term 
(h^/4!)ff  and  are  given  by 

y 


\ ...  1 1  .  -f  ,1  n.  •  n  82  x.  o  sen 
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(3.69) 

i — 1 

II 

i — l 

on 

o 

*5 

C\J 

m 

£3 

+ 

rH 

on 

£3 

+ 

C\J 

o 

£ 

1 — 1 

CM 

£3 

i _ l 

VO 

and 


(3.70) 

6[b2iWi2  +  (b^1+b32)w13]  =  4  . 

The  imposition  of  (3.68)  does  not  reduce  the  number  of 


constraints 

involved  in  a  method  of  rank  three. 

3.3.2 

The  Truncation  Error  The  value  of  t  in 

-  m 

(3.38)  for 

m  =  3  is  obtained  by  finding  the  difference 

between  terms  in  h^/4!  in  (3.62)  and  (3.64).  This  is 


given  by 

(3.71) 

tj  =  u1(D2f+ff  )  +  u2fzDf 

where 


(3.72) 

2  2 

u1  =  1  -  6 (a2wQ2+a^w03) 

and 


(3.73) 

u.  2  ”  1  ™  * 

Following  a  similar  procedure  to  that  used  in 
Section  3.2,  and  employing  (3-47),  the  following  bound 


% 


* 


- 


’ 


V 


' 


can  be  placed  on  the  absolute  value  of 


(3.74)  1 1 3 |  <  (9|u1|+3|u2| )L2M  +  |u1|lm2N_1  . 

It  is  then  possible  to  write  (3.38)  for  RKN  methods  of 
rank  three  as 

(3.75)  | |  <  h^(e1L2M  +  e2LM2N-1)  +  0(h5) 
where 

(3.76)  e1  =  (9 | u1 | +3 | u2 | )/4 ! 
and 

(3.77)  •  e2  =  | u1 | /4 ! 

For  an  ODE  of  type  (3.59),  the  bound  on  the  trun¬ 
cation  error  for  m  =  3  is  given  by 

4 

(3-78)  |Tj  <  7jT(  4  |  un  |  L2M  +  |u.1|LM2N’1)  +  0(h5) 

3  f(x,y)  1  1 

which  was  found  by  again  considering  all  partial  deriva¬ 
tives  with  respect  to  z  as  being  equal  to  zero. 
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CHAPTER  IV 


NEW  SETS  OP  COEFFICIENTS  FOR 
RUNGE-KUTTA-NYSTROM  METHODS  OF  RANK  FOUR 

4 . 1  Two  RKN  Methods 

The  set  of  coefficients  for  RKN  methods  most  widely 
used  for  computational  purposes  is  that  proposed  by 
Nystrom  [16]  and  given  by 


(4.1) 


1/4 

1/4 

1/4 

0 

1/4 

0 

0 

1 

1 

1/2 

1/2 

0 

1/2 

1/2 

0 

0 

1 

1 

1/3 

1/3 

1/3 

0 

1/3 

2/3 

2/3 

1/3 

The  bound  on  the  truncation  error  involved  using  (4.1)  is 

(4.2)  |  T4 1  <  h5(0.0979  L3M  +  0.0354  L2M2N"'1)  +  0(h6) 

calculated  from  (3-56).  The  bound  obtained  for  a  problem 
of  type  (3.59)  is  found  from  (3.61)  and  given  by 


* 
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(4.3)  |T,|  <  h5(0.0111  L3M  +  0.025  L2M2N”1)  +  0(h6) 

4  f(x,y) 

Another  set  of  coefficients  was  suggested  by  Duffin 
[5]  and  given  by 


1/9 

1/9 

1/9 

1/3 

4/9 

1/4 

1/2 

1/4 

1 

1/3 

1/3 

-1/3 

1 

2/3 

1 

-1 

1 

1 

1/4 

1/2 

1/4 

0 

1/4 

3/4 

3/4 

1/4 

where  the  coefficients  e^  and  e^  for  the  truncation 
error  bounds  are 


y" 

el 

e2 

f(x,y,y '  ) 

0.0653 

0.0111 

f (x,y ) 

0.00494 

0.00648 

Application  of  (4.4)  to  several  second  order  ODE's  pro¬ 
duced  better  results  than  those  obtained  using  (4.1)  or 
the  repeated  application  of  any  RK  method  of  rank  four. 


- 
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Superior  stability  was  also  indicated.  These  results 
were  to  be  expected  since  the  error  bounds  for  (4.4)  were 
smaller  than  those  for  (4.1)  [5]. 

4 . 2  The  Two  Sets  of  Constraints 

The  constraints  originally  derived  by  Nystrom  [16] 
and  Duff in  [5]  for  a  RKN  method  of  rank  four  are  given 
by  (3.29),  (3.3D  and  (3.32-37)  which  form  a  system  of 
twenty-seven  non-linear  equations  in  twenty-three  un¬ 
knowns.  As  seven  of  these  equations,  given  by  (3.34), 
(3.36)  and  (3.39),  comprise  the  bound  on  the  truncation 
error,  the  actual  system  necessary  to  be  solved  is  reduced 
to  twenty  equations  in  twenty-three  unknowns.  Attempts 
were  made  to  solve  this  three  parameter  system  but  were 
unsuccessful . 

A  close  examination  of  the  constraints  derived  by 
Nystrom  and  Duffin  revealed  that  if  restriction  (3.30) 
were  imposed,  the  complete  system  would  be  reduced  to 
one  of  twenty-four  equations  in  twenty-three  unknowns 
given  by  (3.29-31).  The  system  to  be  solved  would  not 
include  the  five  constraints  of  (3*39)  and  would  be  a 
four  parameter  system  of  nineteen  equations  in  twenty- 
three  unknowns.  Although  another  degree  of  freedom  was 
obtained,  attempts  to  solve  this  system  were  also 


unsuccessful . 


, 

. 

Of'  .  .  I  19W 


. 


.  .1  3  8  :•  ’>  •  J  a.l 


46 


4 . 3  Two  Procedures  for  Obtaining;  New  RKN  Methods 

Any  set  of  constraints  for  a  RKN  method  of  rank  four 
obtained  from  a  RK  method  of  the  same  rank  using  (2.25) 
satisfies  all  necessary  constraints  derived  by  Nystrom 
and  Duffin.  They  do  not,  however,  satisfy  the  restric¬ 
tions  imposed  by  (3.30)  in  the  new  derivation  of 
constraints . 

4.3.1  Procedure  I  To  obtain  coefficients  for  RKN 
methods  of  rank  four  one  can  obviously  use  (2.25).  Apply¬ 
ing  such  a  RKN  method  to  solve  a  second  order  ODE  is  the 
same  as  using  the  corresponding  RK  method  to  solve  an 
equivalent  system  of  first  order  equations.  To  improve 
upon  these  coefficients,  in  terms  of  truncation  error 
bound  and  the  results  obtained  in  the  solution  of  second 
order  equations,  it  is  noted  that,  in  the  new  derivation, 
only  four  of  the  nineteen  necessary  constraints  contain 
terms  in  b^  .  These  are  given  by 


(4.6) 


a2  b21  » 


(4.7) 


a3  "  b3l  +  b32  > 


(4.8) 


' 


■ 
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and 


(4.9) 


30 a2b  ^ 2w-^ 3  3 0  ( ij 2 ^ 3^  i) 3  ^ ^ 2. 4  —  5  • 


If  one  first  determines  all  coefficients  from  (2.25) 
except  the  b^.  ;  j  <  i,  it  is  possible  to  consider 
(4.7-9)  as  a  two  parameter  system  of  three  equations 

in  the  five  unknowns  b^ia  b325  b4l5  b42  and  b43* 

The  coefficient  b21  can  be  determined  from  (4.6). 

If  b^2  and  t>ij3  are  chosen  as  the  two  parameters 
one  obtains  from  (4.9) 


(4.10) 


1  -  6(a2bi|2+a^bi+3)wll| 


32 


6a2w1^ 


Hence  (4.6-8)  yield 


(4.11) 


2 

a2  a 


(4.12) 


and 

(4.13)  bi,!  =  a4  -  bi)2  _  b43  ‘ 

Using  (4.10-13),  values  for  the  bj,  can  be  obtained 
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which  satisfy  all  necessary  constraints  in  (3.29-31). 

Both  (4.1)  and  (4.4)  can  be  obtained  from  (2.25) 
and  (4.10-13).  Using  the  classical  RK  method  given  in 
Appendix  I,  with  b^  =  0  and  b ^  =  1,  Procedure  I 
gives  (4.1).  Prom  Boulton’s  RK  method  also  in  Appendix 
I,  with  b^  =  and  b^  =  -^-,  the  above  procedure 

yields  (4.4). 

4.3.2  Procedure  II  A  second  procedure  to  determine 
new  RKN  methods  is  to  include  a  fifth  constraint  in  deter¬ 
mining  the  b _  while  again  obtaining  values  for  the 
other  coefficients  as  in  Procedure  I.  The  constraint 

(4.14)  30a2b 32w03  +  30  (a^^+a^b^w^  =  1 

found  in  (3.39)  contains  terms  in  b..  and  the  inclusion 
of  it  enables  (4.7-9)  and  (4.14)  to  be  treated  as  a  one 
parameter  system  of  four  equations  in  five  unknowns. 

Again  b21  can  be  determined  separately  from  (4.6). 

Any  four  unknowns  may  be  expressed  in  terms  of  the 
remaining  one.  If  b^  is  chosen,  (4.9)  and  (4.14) 
yield 

5w0l)  -  wl4 

b32 


(4.15) 


30a2(Wg4Wi3  wQ3w4lt) 


--  ~  M  o-  -4-r 
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The  remaining  b. .  may  be  determined  from  (4.11-13). 

“  J 

The  use  of  either  procedure  reduces  the  bounds  on 
the  truncation  error  given  by  (3.56)  and  (3.61)  compared 
to  that  using  coefficients  obtained  only  from  (2.25). 

The  reduction  is  in  the  term  e^  for  both  bounds  and 
is  larger  using  Procedure  II  since  the  use  of  (4.14)  as 
a  constraint  reduces  | u^ | ,  given  by  (3.43),  to  zero. 

4 . 4  New  Runge-Kutta-Nystrom  Methods 

Consider  the  RK  method  of  rank  four  proposed  by 
Kuntzmann  [3]  and  given  by 


2/5 

2/5 

-3/20 

15/20 

3/5 

19/44 

-15/44 

40/44 

1 

11/72 

25/72 

25/72 

11/72 

Determining  the  corresponding  coefficients  for  a  RKN 
method  from  (2.25)  yields 
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0 

4/25 

3/5 

0 

9/25 

-6/11 

15/11 

0 

1 

2/5 

2/5 

-3/20 

15/20 

3/5 

19/44 

-15/44 

40/44 

1 

11/36 

15/36 

10/36 

0 

11/36 

25/36 

25/36 

11/36 

Any  set  of  coefficients  determined  in  this  manner  will  be 

referred  to  as  the  X  RKN  method  where  X  is  the  name  of 

the  RK  set  of  coefficients;  i.e.  (4.18)  is  the  Kuntzmann 

RKN  method.  If  the  b. .  for  a  RKN  method  are  determined 

ij 

by  either  Procedure  I  or  II,  the  method  will  be  referred 
to  as  a  X  RKN  J  where  X  is  as  defined  previously  and 
J  is  equal  to  I  or  II. 

The  b^j  of  the  Kuntzmann  RKN  I  where  b^2  =  3/5 
and  b^  =  1/5  are  given  by 


4/25 

4/25 

(4.19) 

17/250 

73/250 

90/250 

1/5 

3/5 

1/5 

1 

The  corresponding  b^ 
b^  =  1/5  are 


of  the  Kuntzmann  RKN  II  where 
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4/25 

4/25 

(4.20) 

3/50 

15/50 

18/50 

18/55 

21/55 

11/55 

1 

4.4.1  Truncation  Error  Bounds  The  particular  choices 
for  and  b^  were  made  for  both  Kuntzmann  RKN  I  and 

II  because  they  gave  the  most  substantial  decrease  in  the 
truncation  error  bounds,  not  only  for  Kuntzmann' s  RK  method, 
but  also  for  all  others  used.  This  was  determined  by  try¬ 
ing  different  values  for  the  parameters  within  the  interval 
[0,1]. 

As  all  X  RKN  methods  are  determined  using  (2.25), 

* 

the  truncation  error  bounds  for  these  methods  can  be  con¬ 
sidered  equivalent  to  the  corresponding  bounds  using  the 
X  RK  method.  Because  of  this  one  can  use  this  bound  as 
a  measure  in  comparing  the  error  bounds  of  X  RKN  J 
methods  as  opposed  to  the  use  of  X  RK  methods  in  solving 
a  second  order  ODE. 

The  coefficients  e1  and  e 2  of  the  truncation 
error  bounds  for  Kuntzmann's  RKN  methods,  as  well  as 
Nystrom's,  Duffin' s  and  those  listed  in  Appendix  I,  are 
all  given  in  Table  4.1,  illustrating  the  fact  that,  for 
any  given  X,  e2  is  smallest  for  X  RKN  II  methods. 
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TABLE  4.1 

Truncation  Error  Bounds  for 
Runge-Kutt a-Ny strom  Methods 


Ralston  RKN 
Ralston  RKN  I 
Ralston  RKN  II 
Nystrom  (4.1) 
Classical  RKN 
Classical  RKN  I 
Classical  RKN  II 
Gill  RKN 
Gill  RKN  I 
Gill  RKN  II 
Duffin  (4.4) 
Boulton  RKN 
Boulton  RKN  I 
Boulton  RKN  II 
Kuntzmann  RKN 
Kuntzmann  RKN  I 
Kuntzmann  RKN  II 


f(x,y 

,y’) 

el 

e2 

0.120 

0.0514 

0.120 

0.0169 

0 . 120 

0.0136 

0.0979 

0.0354 

0.0979 

0 . 0458 

0 . 0979 

0.0230 

0.0979 

0 . 0146 

0.0979 

0 .0458 

0.0979 

0.0230 

0 . 0979 

0 . 0146 

0.0653 

0.0111 

0.0653 

0.0417 

0.0653 

0 . 0111 

0.0653 

0.00694 

0.055 

0.0333 

0.055 

0 . 018 

0.055 

0 . 01 

f(x,y) 


e2 


0.0121 

0.0287 

0.0121 

0.0112 

0.0121 

0.00907 

0.0111 

0.025 

0 . 0111 

0.025 

0.0111 

0.0167 

0 . 0111 

0.00833 

0.0111 

0.025 

0.0111 

0.0167 

0 . 0111 

0.00833 

0.00494 

0.00648 

0.00494 

0.0222 

0 . 00494 

0 . 00648 

0.00494 

0.00370 

0 . 00889 

0.0167 

0 . 00889 

0.012 

0.00889 

0.00667 

' 
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A  close  examination  of  Table  4.1  reveals  there  is 
no  method  that  always  has  the  lowest  values  of  both  e^ 
and  e2.  Using  the  sum  of  e^  and  e^  as  a  criterion, 
although  this  is  not  necessarily  significant,  it  is 
obvious  that  among  the  methods  listed  in  Table  4.1  the 
Kuntzmann  RKN  II  method  has  the  smallest  bound  when 
y"  =  f(x,y,y’),  and  the  Boulton  RKN  II  method  for  f(x,y). 


. 


CHAPTER  V 


NUMERICAL  RESULTS 


5 • 1  Test  Problems 

To  assess  the  relative  merits  and  demerits  of  the 
X  RKN  J  methods  in  solving  second  order  ODE's  numerical 
solutions  of  several  problems  with  known  analytical  solu¬ 
tions  were  computed  by  use  of  both  RK  and  RKN  methods. 
These  test  problems  are 


(5.1) 


y"  =  y;  y(0)  =  y’(0)  =  l 


(5.2) 


y"  -  -0.5  +  1.5y' ;  y(0)  =  y' (0)  =  1 


(5.3) 


y"  =  lly  -  lOy ' ;  y(0)  =  y'(0)  =  1  , 


(5. 4) 


y"  =  y;  y(o)  =  l,  y'(o)  =  -l 


(5.5) 


y"  =  -0.5y  -  l.5y';  y(0)  =  l,  y'(0)  =  -l  , 


(5.6) 


y"  =  lly  +  10y';  y(0)  =  1,  y'(0)  =  -1 


(5.7) 


y"  =  y;  y(0)  =  2,  y'(0)  =  0  , 


(5.8) 


y 


n 


y  -  2  sin  x;  y(0)  =  1,  y’(0)  =  0  , 


' 


■ 
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^5.9)  y"  =  -y;  y(0)  =  0,  y'(0)  =  1 


and 


(5.10)  y"  =  -y;  y(0)  =  1,  yf(0)  =  0  . 

Their  corresponding  analytical  solutions  are  given  by 
y  =  ex  for  (5.1-3);  y  =  e~x  for  (5.4-6);  y  =  ex  +  e“x 
for  (5-7);  y  =  sin  x  +  e”x  for  (5.8);  y  =  sin  x  for 
(5.9);  and  y  =  cos  x  for  (5.10).  The  above  problems 
are  often  used  in  the  literature  [3],  [4]  and  most  have 
exponential  solutions  since  this  type  of  problem  is  best 
suited  to  measure  accuracy  and  stability. 

All  problems  are  solved,  both  as  a  simultaneous 
system  of  first  order  equations  using  Ralston’s  RK  method 
of  rank  four,  and  as  a  second  order  ODE  by  Ralston  RKN  II 
and  Kuntzmann  or  Boulton  RKN  II  methods  of  corresponding 
rank  (Appendix  I).  These  particular  methods  are  used 
because  Ralston  RK  (Rl)  method  has  the  minimum  truncation 
error  bound  for  a  first  order  equation  while  Kuntzmann 
RKN  II  (K2)  has  the  corresponding  bound  for  second  order 
equations  where  yM  =  f(x,yay’),  and  Boulton  RKN  II  (B2) 
for  yM  =  f(x,y).  The  numerical  results  are  listed  in 
Appendix  III. 
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Comparisons  in  accuracy  and  stability  are  made  for 
all  test  problems  by  listing  the  ratios  of  the  actual 
errors  involved  for  each  method.  In  this  manner  the  RKN 
methods  are  compared  with  the  RK  methods  used  as  well  as 
with  each  other. 

Several  other  initial  value  problems  with  analytical 
solutions  involving  logarithmic  and  trigonometric  functions 
were  solved.  The  results,  although  not  listed  in  Appendix 
III,  favoured  the  use  of  the  X  RKN  J  methods. 

5 . 2  Accuracy 

The  results  of  problems  (5.1),  (5.4)  and  (5.7-10) 
indicate  the  superiority  of  RKN  methods  over  RK  methods 
in  terms  of  initial  accuracy.  For  these  six  problems 
the  X  RKN  J  methods  give  up  to  three  more  significant 
digits  of  accuracy,  than  the  Ralston  RK  method;  i.e. 
problem  (5.1),  (5.4),  (5.7)  and  (5.10). 

The  Boulton  RKN  II  method,  used  for  the  above 
mentioned  problems  because  all  are  of  the  form 
y"  =  f(x,y),  gives  better  results  than  the  Ralston 
RKN  II  method.  For  (5.7)  and  (5.10)  the  difference  is 
as  much  as  two  significant  figures  initially  while  for 
(5.1),  (5.4)  and  (5.8-9)  the  difference  is  not  as  signifi¬ 
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The  second  order  ODE's  given  by  (5.2-3)  and 
(5.5-6)  are  of  the  form  y"  =  f(x,y,y'),  hence  the 
Kuntzmann  RKN  II  method  was  used.  These  problems  were 
specifically  constructed  to  illustrate  what  can  occur 
in  solving  differential  equations  numerically.  In  all 
four  problems  the  Ralston  RK  method  gives  better  results 
than  the  X  RKN  II  methods  used.  A  close  examination  of 
the  computational  procedure  followed  for  these  problems 
indicates  the  results  can  be  attributed  to  a  significant 
build-up  in  round-off  error.  This  is  caused  by  the  sub¬ 
traction  of  two  numbers  of  approximately  the  same  magni¬ 
tude  and  occurs  for  the  RKN  methods  because  of  the 
restrictions  on  the  coefficients  imposed  by  (3.30). 

The  propagation  of  round-off  error  is  magnified  in  these 
four  problems  because  of  the  form  of  the  function 
fCxjyjy*)*  and  is  especially  noticeable  in  the  solu¬ 
tions  of  (5*3)  and  (5.6). 

The  Ralston  RKN  II  method  gives  better  results  for 
(5.2-3)  while  Kuntzmann' s  is  superior  for  (5.5-6).  How¬ 
ever,  the  differences  in  these  results  are  insignificant. 

5 . 3  Stability 

Problems  (5.1)  and  (5.4)  were  integrated  over  a  larger 
interval  than  other  problems  to  indicate  the  stability  of 
the  different  methods.  The  results  for  (5.1)  indicate  not 
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only  superior  accuracy  for  the  X  RKN  II  methods,  but 
this  accuracy  is  maintained  throughout  the  interval 
of  integration.  Conversely,  although  the  solution  for 
(5.4)  is  initially  more  accurate,  the  accumulation  of 
round-off  error  is  significant  and  results  in  instability 
for  this  problem  using  X  RKN  II  methods.  The  Ralston  RK 
method  maintains  stability  throughout  the  interval  as  in 
(5.1). 

The  significance  of  the  round-off  error  in  (5.4) 
can  be  attributed  to  restriction  (3*30)  and  the  fact  that 
y,  y’  and  y"  are  of  equal  absolute  value.  This  causes 
the  subtraction  of  numbers  of  approximately  the  same 
magnitude . 

The  instability  as  well  as  inaccuracy  for  the  solu¬ 
tions  of  (5*2-3)  and  (5.5-6)  can  also  be  attributed  to 
restriction  (3.30),  and  the  values  of  y,  y’  and  y". 

The  form  of  the  function  f(x,y,yf)  magnifies  significantly 
the  accumulation  of  error  in  these  problems.  The  obvious 
differences  in  accuracy  and  stability  for  (5.1-3)  and 
(5*4-6)  illustrate  the  sensitivity  of  the  solution  to 
the  form  of  the  function  f(x,y,y').  These  four  problems 
maintain  stability  using  the  Ralston  RK  method. 

The  ODE's  given  by  (5.7-10)  have  stable  as  well  as 
accurate  solutions  throughout  the  interval  of  integration 
for  both  the  Ralston  RK  and  the  X  RKN  II  methods  used. 


. 


'■ 

, 

’ 


■ 


59 


There  was  no  particular  distinction  between  the  stability 
of  the  X  RKN  II  methods  used  for  the  test  problems. 


« 
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CHAPTER  VI 


CONCLUSIONS  AND  SUGGESTIONS 
FOR  FURTHER  RESEARCH 


6 . 1  Conclusions 

The  theoretical  and  numerical  results  of  this  study 
lead  to  the  following  conclusions : 

1.  The  X  RKN  J  methods  of  rank  four  obtained  from 
Procedures  I  and  II  are  superior  to  the  use  of 
the  corresponding  RK  method  of  the  same  rank 
with  regards  to  the  truncation  error  bound  found 
by  the  approach  of  Lotkin  [15]. 

2.  The  X  RKN  II  methods  are  generally  more  accurate 
than  the  corresponding  RK  methods  of  the  same 
rank . 

3.  The  X  RKN  J  methods  appear  to  be  more  susceptible 
to  the  propagation  of  round-off  error  than  the 

RK  methods.  In  the  sample  problems  tested  this 
was  attributed  to  the  added  restrictions  imposed 
on  the  coefficients  by  (3.30)  and  the  fact  that 
the  absolute  magnitudes  of  ys  yT  and  yM  were 
equal . 

4.  The  use  of  either  a  RK  or  RKN  method  of  cor¬ 
responding  rank  involves  the  same  number  of 
operations  and  function  evaluations  for  each 
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integration  step  if  the  methods  are  programmed 
efficiently  with  respect  to  time  and  round-off 
error . 

5.  No  single  X  RKN  II  method  appears  best  suited 
for  all  problems  although  the  best  results  were 
obtained  from  the  methods  attributed  to  Ralston, 
Kuntzmann  and  Boulton.  The  method  attributed  to 
Boulton  is  that  best  suited  for  functions  of  the 
form  f(x,y). 

6.  RKN  methods  have  all  the  advantages  and  disadvant¬ 
ages  listed  in  Section  1.3  for  RK  methods  and  thus 
are  best  used  only  in  obtaining  starting  values 
for  the  solution  by  predictor-corrector  methods. 

6 . 2  Suggestions  for  Further  Research 

1.  The  general  relationships  between  coefficients 
of  RK  and  RKN  methods  of  the  same  rank  given  by 
(2.25)  provide  the  basis  for  continuing  the 
research  previously  described  for  RKN  methods 
of  degree  greater  than  four.  Perhaps  a  general 
procedure  could  be  developed  to  obtain  modified 
coefficients  for  a  RKN  method  of  any  degree. 

2.  The  general  direct  Runge-Kutta  type  methods  for 
third  and  higher  ordered  ODE’s  could  be  derived 

in  the  same  manner  as  the  RKN  method  in  Chapter  II. 
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Relationships  could  then  be  established  between 
the  coefficients  for  these  methods  and  for  the 
RK  methods  for  first  order  equations. 

3.  A  third  or  higher  ordered  equation  can  be  treated 
as  any  system  of  lower  order  ODE’s,  hence,  direct 
RK  type  methods  could  not  only  be  derived  from  a 
system  of  first  order  equations  as  in  Chapter  II, 
but  also  from  some  system  of  ODE’s  containing 
equations  of  higher  order  than  one.  For  example, 
a  third  order  equation  could  be  treated  as  a 
system  consisting  of  a  first  and  second  order 
ODE  while  a  fourth  order  equation  could,  perhaps, 
be  treated  as  a  pair  of  second  order,  or  first 
and  third  order  ODE’s.  In  this  manner  general 
relationships  could  be  derived  between  the  coef¬ 
ficients  ,of  different  RK  type  methods. 

4.  The  two  parameter  non-linear  system  for  RK 

methods  of  rank  four  has  been  solved,  as  mentioned 
in  Chapter  I,  in  terms  of  and  Hence, 

Procedure’s  I  and  II  have  effectively  solved  the 
non-linear  system  necessary  for  coefficients  of 
RKN  methods  of  rank  four  in  terms  of  four  and 
three  parameters  respectively.  Perhaps  one  could 
express  the  coefficients  of  RKN  methods  of  rank 
four,  as  well  as  methods  for  higher  ordered 
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equations  of  the  same  rank,  in  terms  of  a  few 
parameters.  This  would  obviously  enable  a  more 
comprehensive  study  of  such  methods. 
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APPENDIX  I 


SOME  RUNGE-KUTTA-NYSTROM  METHODS 


OP  RANK  POUR 


This  appendix  contains  the  coefficients  for  the  best 

known  and  most  widely  used  RK  methods  of  rank  four  for 

first  order  ODE’s,  and  the  corresponding  RKN  coefficients 

obtained  using  (2.25),  Procedure  I  and  Procedure  II. 

The  methods  listed  are  the  Classical  and  those  due 

to  Boulton,  Gill  [10]  and  Ralston  [ 1 8 ] .  The  RK  methods 

are  presented  first  and  labelled  RK  while  the  coefficients 

obtained  using  only  (2.25)  are  listed  under  RKN.  The 

b.  found  from  Procedures  I  and  II  are  labelled  PI  and 

PII,  respectively,  and  listed  along  with  the  corresponding 

2 

values  of  the  a^ . 

The  RKN  methods  presented  by  Nystrom  [16]  and  Duffin 
[5]  are  given  in  Chapter  IV  with  the  RK  and  RKN  methods 
attributed  to  Kuntzmann  [3]. 


' 


■ 

■ 


■  I 


69 


CLASSICAL 


RK 


0  .  5 

0.5 

0 

0  .  5 

0.  5 

0 

0 

1 

1 

0.16666667 

0. 33333333 

0.33333333 

0. 16666667 

RKN 


0 

0.25 

0.  5 

0 

0.25 

0 

1 

0 

1 

0.  5 

0.5 

0 

0 . 5 

0.5 

0 

0 

1 

1 

0. 33333333 

0.33333333 

0. 33333333 

0 

0.33333333 

0. 66666667 

0.66666667 

0. 33333333 

PIT 


hi  ■  r  kM 
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BOULTON 


RK 


0. 33333333 

0. 33333333 

"0. 33333333 

1 

0.66666667 

1 

“1 

1 

1 

0.125 

0.375 

0.375 

0.125 

RKN 


0 

0.11111111 

0.66666667  0 

0.44444444 

“1.3333333  2  0 

1 

0. 33333333 

0. 33333333 

“0.33333333  1 

0. 66666667 

1  “1  1 

1 

0.25  0.5  0.25  0 

0.25  0.75  0.75  0.25 


PI 


0. 11111111 

0.11111111 

0 .11111111 

0. 33333333 

0 . 44444444 

o 

• 

ro 

0. 6 

CM 

• 

O 

1 

PII 


0.  11111111 

0.11111111 

0 . 044444444 

0.4 

0.44444444 

0.4 

0.4 

CM 

• 

o 

1 

.0 


V 
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GILL 


RK 


0 . 5 

0 . 5 

0 . 20710678 

0 .29289322 

0 

0. 5 

0 

"0.70710678 

1. 7071068 

1 

0.  16666667 

0. 097631073 

0 .56903559 

0. 16666667 

RKN 


0 

0.25 

0.29289322 

0 

0.25 

0 

1 

0 

1 

0.  5 

0.5 

0 .20710678 

0.29289322 

0.5 

0 

"0. 70710678 

1 .7071068 

1 

0 . 33333333 

0.097631073 

0.56903559 

0 

0. 33333333 

6. 19526215 

1 .1380712 

0.33333333 

PI 


0.25 

0.25 

0. 19142136 

0 . 058578644 

0.25 

o 

• 

ro 

o 

• 

CD 

0.2 

1 

PIT 


0 .25 

0.25 

0.13284271 

0.11715729 

0.25 

0.4 

0.4 

0.2 

1 
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RALSTON 


RK 


0.4 

0.4 

0.2969776 

0. 15875966 

0 .45573726 

0. 21810038 

"3.0509647 

3. 8328643 

1 

0 . 17475028 

"0. 55148053 

1. 2055355 

0.17118478 

RKN 


0 

0.16 

0.  12700773 

0 

0 .20769645 

"0. 16422207 

1 . 2170085 

0 

1 

0.4 

0.4 

0 . 2969776 

0. 15875966 

0.45573726 

0. 21810038 

"3. 0509647 

3. 8328643 

1 

0. 34952057 

"0. 66177664 

1 .3122561 

0 

0 . 34952056 

"1. 1029611 

2.4110709 

0. 34236956 

PI 


0 .16 

0.  16 

0. 15243899 

0.055257457 

0.20769645 

0.2 

0.6 

0.2 

1 

PIT 


0 .16 

0.16 

0. 14419258 

0.  063503866 

0. 20769645 

0.25807373 

0. 54192627 

0.2 

1 
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APPENDIX  II 


APL  FUNCTION  LISTINGS 

APL  is  a  machine  executable  version,  available  in 
a  time-sharing  environment,  of  what  is  more  commonly 
known  as  the  Iverson  Language.  An  interpreter  has  been 
written  for  the  IBM  360,  Models  50  and  higher.  This 
particular  language  was  chosen  for  the  numerical 
experiments  because:  1)  all  computations  are  performed 
in  double  precision;  2)  the  time-sharing  environment 
provided  a  means  of  comprehensively  testing  all  methods 
and  problems,  comparing  results  as  they  were  obtained. 
For  a  complete  introduction  to  this  language  see 
Introduction  to  APL  360  Concepts  by  W. S.  Adams,  Depart¬ 
ment  of  Computing  Science,  University  of  Alberta. 

The  following  'pages  contain  listings  of  the  more 
important  functions  used  in  testing  the  X  RK  and 
X  RKN  J  methods. 

USE  1  X  and  USE2  X  convert  the  coefficients  for 
RK  and  RKN  methods,  respectively,  into  the  proper  vector 
and  matrix  form  necessary  for  use  in  solving  a  problem. 
The  parameter  X  is  a  vector  of  length  thirteen  for 
USE  1  and  of  length  twenty-three  for  USE 2,  the  number 
of  coefficients  not  equal  to  zero  in  RK  and  RKN  methods, 
respectively,  of  rank  four. 
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X  RK 1  Y  will  solve  any  simultaneous  system  of 
first  order  equations  using  a  RK  method  of  any  rank. 
YIPRIME  is  used  to  evaluate  the  functions  of  the 
system  being  solved  and  YEXACT  to  calculate  the 
analytical  solution. 


X  RK  2  y 

solves  second  order  ODE’s  directly  using 

RKN  methods . 

Y2PRIME  is  used  to  evaluate  the  function 

involved  and 

YEXACT  to  calculate  the  analytical  solution. 

RATIO  X 

calculates  the  ratio  of  the  actual  errors 

involved  in  the  different  methods.  These  errors  have 
been  previously  calculated  and  saved  in  RK 1  and  RK 2. 

The  parameter  X  used  in  RK 1 ,  RK 2  and  RATIO  is 
a  four  element  vector  where  ZCID^Xq,  y[2]«-h,  X[3]«- 
Man  upper  limit  on  the  interval  of  integration",  and 
y[4]+-m,  the  rank  of  the  method. 

The  parameter,  Y  is  a  vector  of  length  n  in 
RK 1  where  n  is  the  number  of  equations  in  the  system 
being  solved,  and  of  length  two  in  RK 2,  i.e.  y+y0,  ^0* 


. 
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vusEiiniv 

V  USE  1  X 

[1]  B1+-  4  4  pAl«-J/l«-4p  0 

[2]  i41[l+i3>*[  i  3] 

C3]  S1[2;1>X[4] 

[4]  25 1  [  3  ;  i2>*[4  +  i2] 

[5]  5 1  [  4  ;  i  3  ]  +-X  [  6  +  \  3  ] 

[6]  Wll\  4]«-Ar[  9  +  i  4] 

V 


lUSE2imv 

V  USE  2  X 

[1]  B  2*-C2*-  4  4  p  W2-*-  2  4  p42«-4pO 

[2]  A2Z  l+i  3>*[  i  3] 

[3]  5  2  [  2  ;  1  ]«-J[  4  ] 

[4]  5  2  [  3  ;  i  2  ]  <-XL  4  +  i  2  ] 

[5]  B2[4  ;  i  3]«-*[6+i  3] 

[6]  C2[2  ;  1>X[  10] 

[7]  C  2  [  3  ;  i  2  ]  *-X  [  1  0  +  i  2  ] 

C8]  C2[4;  i  3]«-X[12  +  i  3] 

[9]  W2L  x  2  ;  ]<-  2  4  p*[15+i8] 

V 


VBXIED] V 

V  M+X  RK 1  XiIiJ-,K;L 

[1]  *[3>U[3]-Z[i])i*[2] 

[2]  SAVE+M+\I+0 

[3]  K+i  (  p  Y )  ,  XL  )p  J<-0 

[4]  Kl  ;t7]^-Z[2]x(  J[l]+i4l[J-]xJ[2]  )  Y1PRIME  (  J+  +  /  (  (  (pY)  tXL 
4]  )  p  B 1  [  («7  +-J  +1)  ;  ]  )*K) 

[5]  -U<*[ 4])/4 

[6]  I«-I+l 

[7]  M+M,XH1 ,7,  |  (  (pY)YEXACT  XL  1  ]«-*[  1 1  +  XL  2  ]  ) -Y+Y+  +  /  (  (  (  p  7)  , 
*[4] )pWl)xK 

[8]  +  U[3]>I)/3 

[9]  M«-U[3]  ,  l  +  2xp7)pM 

[10]  ML  ;  3  5>*[2]xAf[;  3  5] 

[11]  S>17£*-S>ira,Af[  S4]  ,ML  ;5] 

[12]  (6p»  ')t'X',(12p'  »),'J,,(12  p»  *  )  •  *  ^  *  .  (  10o  1  ,),’|£'(ir) 

I  '  , ( 7p  »  •  )  ,  1  \E(  V)  |  ' 


V 
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Vi?Z2[D]  V 

V  M+X  RK2  Y;I;J;K;LiXViYV 

[1]  *[3>U[3]-*[1]  )*XZ2l 

[2]  SAVE+M+\ 1+0 

[3]  K-+-XL  4  ]  pt7-<-0 

[4]  XV+XH1+A21J+-J+13XX121 

[5]  YV+(YZ  1]  +  U2[c7]x*[2]xy[2]  )  +  +  /B  2  [  J  ;  ]  *K  )  ,  (  Yl  2  ]  +  +  / ( 
2t*[2] )xC2[J; ]x*) 

[6]  KlJ]+((Xl2]*2)i2)*XV  Y2PRIME  YV 

[7]  ->(J<Z[4]  )/4 

[8]  y[l>y[l]  +  U[2]xJ[2]  )  +  +  /V2[l;]xZ 

[9]  y[2]«-y[2]+  +  /0/2[2;]xJOm2] 

[10]  J«-I+l 

[11]  YV+(pY)YEXACT  XL  1  ]«-*[  1  ]  +  *[  2  ] 

[12]  M«-M,*[i],y[i]tU[2]xy[2]),(  imi]-i[i]),IU[ 

2 ] x  J7[ 2 ] )-J[2]xJ[2] 

[13]  -U[3]>I)/3 

[14]  M+(XL3l95)pM 

[15]  SAVE+SAVE  ,ML  ;  4]  ,Ml  ;  5] 

[is]  ( 6 p *  •),,y’>(i2P'  •),,y,,(i2p'  M.'F’.ciop'  'K'ltfcn 
I  '  ,  (  7p  '  »  )  ,  '  |£(  V)  |  » 

V 


VRATIOl Q]V 


[1] 

[2] 

[3] 

[4] 

[5] 

[6] 

[7] 

[8] 

[9] 

[10] 
[11] 
[12] 

[13] 

[14] 

[15] 


M+RATIO  X;I 
*[4>U[3]-Z[1]  )*XZ  2] 

Sl^SrSl 
S2+SiS2 
S+S2iSl 
M+(XZ 4] , 6  ) p  0 
M[;l>Sl[iZ[4]] 

MZ  ;2>Sl[Z[4]  +  iZ[4]] 

MZ  ;3>S2[ iZ[4]] 

MZ  ;4]<-S2[Z[4]  +  iX[4]] 

MZ  ;5>S[i*[4]] 

MZ  ;6>S[Z[4]  +  iZ[4]] 

(  2 6p '  •  ) , 'RATIO  OF  ERRORS' 

I  r 

(  6  p  ’  '  )  ,'E(R1)/E(B2)'  AHP'  ' 
lip'  ' ) , 'E(B2)/E(R2 ) ' 

( 5  p  ’  '),'7',(10p'  ’),’7'f(10p 

lOp ’  •  ) ,  'Y'  ,(10p '  '  ) , ’ V' 


E(R1) /E(R2) ' ,( 

) , '  y ' ,  ( i op  *  » ) , '  V' ,  ( 


v 
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APPENDIX  III 


LISTING  OP  NUMERICAL  RESULTS 

This  appendix  contains  the  numerical  results  found 
using  the  Ralston  RK  (Rl)  and  RKN  II  (R2)  methods  to 
solve  a  second  order  ODE.  The  Boulton  (B2)  or  Kuntz- 
mann  (K2)  RKN  II  method  is  also  used  in  solving  the 
problem.  The  results  are  compared  by  means  of  a  table 
of  the  ratio  of  the  actual  errors.  Listed  are  the  ratios 
of  the  errors  involved  in  the  Ralston  RK  method  over  the 
Ralston  RKN  II  and  the  Boulton  or  Kuntzmann  RKN  II 
methods;  and  the  Boulton  or  Kuntzmann  RKN  II  method 
over  the  Ralston  RKN  II  method. 

The  APL  functions  Y1PRIME,  Y2PRIME  and  YEXACT  are 
listed  with  the  results,  as  well  as  the  parameter  X 
and  Y  of  the  functions  RK 1,  RK 2  and  RATIO  (Appendix 
II).  The  problems  are  solved  in  the  order  of  listing  in 
Chapter  V;  i.e.  (5.1-10).  Most  are  solved  over  the 
interval  [0,2]  although  (5-1)  and  (5.*0  are  integrated 
from  [0,6].  The  interval  h  is  in  all  cases  equal  to 


0.2. 


1 


' 


! )  t 

• :  .  j 


. 


78 


(5.1) 

vyippiM#[[]]v 

V  V+X  yiPPJM#  X 

[1]  F«*y[  2]  ,  J[  1] 

V 

X 

0  0.2  6  4 

J 

1  1 


USE  1  PI 
X  RK 1  J 


X 

y 

V 

1 E(X) | 

1  #(  F)  | 

2. 00000#“l 

1 . 22140P0 

2. 44280P~1 

2  .  75818#“ 

6 

S.  51635# 

4. 00000#“l 

1 . 49182P0 

2. 9  8  3  6  4P~ 1 

6  .  73768#“ 

6 

1 .34754# 

6 . 00000#“l 

1 . 82211P0 

3 . 64421P~1 

1 .23441#“ 

5 

2 . 46882# 

8 . 0  0  0  0  OP  1 

2 . 22552P0 

4.4510  4P  1 

2. 01028#“ 

5 

4.02057# 

1 . OOOOOPO 

2 . 71825P0 

5 . 4  3  6  5  0P~1 

3 .06920#“ 

5 

6.13841# 

1  .  2  0  0  0  OP  0 

3 . 32007P0 

6 . 64014P"1 

4 .49848#“ 

5 

8.99695# 

1 . 40000P0 

4. 05514P0 

8 . 1102 1E~ 1 

6 .41018#“ 

5 

1 .28204# 

1 . 6  0  0  0  OP  0 

4. 95294P0 

9 . 90589P"l 

8 .94789#“ 

5 

1.78958# 

1.  8  0  0  0  0#0 

6 . 04952P0 

1.  2  0  9  9  OP  0 

1 .22951#“ 

4 

2.45902# 

2. OOOOOPO 

7 . 38889P0 

1. 47778P0 

1 . 66858#“ 

4 

3.33717# 

2  .  2  0  0  0  OP  0 

9 . 02479P0 

1 . 80496P0 

2 .24181#“ 

4 

4.48362# 

2 . 4  0  0  0  OP  0 

1 . 10229P1 

2. 20458P0 

2 . 98707#“ 

4 

5 .97414# 

2.  6  0  0  0  OP  0 

1 . 34633P1 

2  .  6  9  2  6  IE  0 

3 . 95245#“ 

4 

7.90490# 

2.  8  0  0  0  OP  0 

1 . 64441P1 

3 . 28883P0 

5 . 19887#“ 

4 

1 . 03977# 

3  .  OOOOOPO 

2 . 00849P1 

4 . 01697P0 

6 . 80348# 

4 

1.36070# 

3 . 2  0  0  0  OP  0 

2 . 45316P1 

4. 90633P0 

8 .86376#“ 

4 

1 . 77275# 

3 . 40000P0 

2 . 99629P1 

5 . 99259P0 

1 .15029#“ 

3 

2 .30057# 

3 . 6  0  0  0  OP  0 

3 . 65967P1 

7 . 31935P0 

1 . 48760#“ 

3 

2  .  97521# 

3 . 80000P0 

4 . 46993#1 

8 . 93985P0 

1 .91790#“ 

3 

3 . 83581# 

4. OOOOOPO 

5 .45957#! 

1. 09191P1 

2 .46582# 

3 

4.93164# 

4 . 2  0  0  0  OP  0 

6 . 66832#! 

1 . 33366P1 

3 .16235#“ 

3 

6. 32469# 

4 . 4  0  0  0  OP  0 

8 . 14468P1 

1 . 62894P1 

4 .04642#“ 

3 

8 .09285# 

4 . 6  0  0  0  OP  0 

9 . 94791P1 

1 . 98958P1 

5 . 16696#“ 

3 

1 . 03339# 

4 . 8  0  0  0  0#0 

1 . 21504P2 

2 . 43  008P1 

6 . 58532#“ 

3 

1 .31706# 

5 . OOOOOPO 

1 . 48405P2 

2 . 9  6  81  OP 1 

8  .  37845# 

3 

1 .67569# 

5  .  2  0  0  0  OP  0 

1 . 81262P2 

3 . 62523P1 

1 .06428#“ 

2 

2. 12856# 

5 . 4  0  0  0  OP  0 

2 . 21393P2 

4 . 42786P1 

1 .34991#“ 

2 

2.69982# 

5  .  6  0  0  0  OP  0 

2 . 70409P2 

5 . 40819P1 

1 .70985#“ 

2 

3 .41969# 

5  .  8  0  0  0  OP  0 

3 . 30278P2 

6 . 60556P1 

2.16299#“ 

2 

4. 32599# 

6 . OOOOOPO 

4 . 03401P2 

8 . 06803P1 

2 .73298#“ 

2 

5 . 46597# 

7 

6 

6 

6 

6 

6 

5 

5 

5 

5 

5 

5 

5 

4 

4 

4 

4 

4 

4 

4 

4 

4 

3 

3 

3 

3 

3 

3 

3 

3 
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(5.1) 

1Y2PRIMEIU ]V 

V  V+X  Y2PRIME  Y 

Cl]  y«-y[l] 

V 


USE  2  5  2 

y  552  y 


X 

y 

V 

I5(y) | 

|5(7)| 

2  .  000005"l 

1 . 2214050 

2 . 442815~1 

2 .604615" 

9 

7 .276485" 

4.0000  0E~ 1 

1. 4918250 

2. 983655~1 

1 .316325" 

8 

1 . 641495" 

6.000005  1 

1. 8221250 

3.644245  1 

3 .383765" 

8 

2.812965" 

8 .000005"! 

2. 2255450 

4.451085  1 

6. 757995" 

8 

4.331525" 

1 . 0000050 

2. 7182850 

5.436565  1 

1 . 183375" 

7 

6. 309975" 

1.2000050 

3. 3201250 

6 . 640235"1 

1 .913135" 

7 

8. 891065" 

1 .4000050 

4. 0552050 

8.110405  1 

2. 933045" 

7 

1.22557 5" 

1 . 6000050 

4. 9530350 

9 . 906065"1 

4. 331285" 

7 

1 .663355" 

1 . 8000050 

6 . 0496550 

1 . 2099350 

6. 221665" 

7 

2.231545" 

2. 0000050 

7 . 3890650 

1 .4778150 

8.75052 5" 

7 

2. 967055" 

2. 2000050 

9 . 0250150 

1 . 8050050 

1 . 210545" 

6 

3 .916615" 

2 .4000050 

1 . 1023251 

2 . 2046350 

1 .652625" 

6 

5 .139365" 

2. 6000050 

1 . 3463751 

2 .6927550 

2. 231865" 

6 

6.71005 5" 

2. 8000050 

1 . 6444651 

3 .2889350 

2 . 987205" 

6 

8 .723095" 

3 . 0000050 

2. 0085551 

4.0171150 

3 .968125" 

6 

1 .129765" 

3 . 2000050 

2. 4532551 

4. 9065050 

5 .237365" 

6 

1.458375" 

3 .4000050 

2. 9964151 

5 . 9928250 

6 . 874475" 

6 

1 .877045" 

3 . 6000050 

3 . 6598251 

7 . 3196450 

8 . 980035" 

6 

2 . 409555" 

3 . 8000050 

4. 4701251 

8 .9402350 

1 . 168125" 

5 

3 .085815" 

4 .0000050 

5 .4598151 

1 .0919651 

1. 513855" 

5 

3.943375" 

4. 2000050 

6 . 6686351 

1 . 3337351 

1 . 955425" 

5 

5 .029365" 

4.4000050 

8 . 1450851 

1 .6290251 

2 . 518325" 

5 

6 .402925" 

4 .6000050 

9 . 9484351 

1 .9896951 

3. 234645" 

5 

8.138145" 

4. 8000050 

1. 2151052 

2 .4302151 

4.144715" 

5 

1 .032785" 

5 . 0000050 

1. 4841352 

2.9682651 

5 .299215" 

5 

1 .308825" 

5 . 2000050 

1 . 8127252 

3 . 6254451 

6 . 761775" 

5 

1 . 656465" 

5 . 4000050 

2. 2140652 

4 .4281351 

8. 612205" 

5 

2.093885" 

5 . 6000050 

2 . 7042652 

5 .4085351 

1 . 095065 

4 

2. 643795" 

5 . 8000050 

3 . 3029952 

6 .6059951 

1.  390225" 

4 

3.334575" 

6 . 0000050 

4. 0342952 

8 .0685751 

1 .  762405" 

4 

4. 201635" 

9 

8 

8 

8 

8 

8 

7 

7 

7 

7 

7 

7 

7 

7 

6 

6 

6 

6 

6 

6 

6 

6 

6 

5 

5 

5 

5 

5 

5 

5 


S1<-SA  VE 


' 
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(5.1) 

USE  2  R 2 
X  RK  2  Y 


X 

7 

V 

I2(y) | 

|2(7)| 

2 . 0  0  0  0  02”l 

1 . 2214020 

2. 4428 12”l 

1 .515822” 

8 

1 . 490632 

4 . 0  0  0  0  0E  1 

1 .4918220 

2.98365 2  1 

3 . 247602” 

8 

2.73082 2 

6 .  0  0  0  0  0E~ 1 

1 . 8221220 

3.644242  1 

5 .299202” 

8 

3 .701662 

8 . 0  0  0  0  0E~ 1 

2. 2255420 

4.451082  1 

7 . 794902” 

8 

4. 358202 

1 . 0000020 

2 .7182820 

5 . 436562”l 

1 .088612" 

7 

4. 624322 

1  .  2  0  0  0  OE 0 

3 . 3201220 

6 . 640232“l 

1 .475952“ 

7 

4.385532 

1 . 4  0  0  0  OE 0 

4. 0552020 

8.110402  1 

1 .964682” 

7 

3 .479372 

1 . 6  0  0  0  0E0 

4 .9530320 

9.906062  1 

2 . 583782“ 

7 

1 .682752 

1 . 8  0  0  0  OE 0 

6 .0496520 

1. 2099320 

3 . 369472” 

7 

1 .304502 

2 . 0000020 

7 . 3890620 

1.4778120 

4. 367242” 

7 

5 . 880842 

2 . 2  0  0  0  OE 0 

9 .0250120 

1. 8050020 

5. 634122” 

7 

1 .257022 

2. 4  0  0  0  02?  0 

1 . 1023221 

2. 2046420 

7. 241752” 

7 

2. 205652 

2 . 6000020 

1 . 3463721 

2. 6927520 

9 . 280022“ 

7 

3.52282 2 

2  .  8  0  0  0  02? 0 

1 . 6444621 

3. 2889320 

1 . 186172 

6 

5 . 323372 

3 . 0000020 

2 . 0085521 

4. 0171120 

1 . 512842“ 

6 

7.755302 

3.2000  0E0 

2.4532521 

4. 9065120 

1 .925732“ 

6 

1. 100872 

3. 4  0  0  0  OE 0 

2 .9964121 

5. 9928220 

2.447062“ 

6 

1. 532712 

3 . 6000020 

3 . 6598221 

7 . 3196520 

3 . 104662” 

6 

2.102172 

3. 8000020 

4. 4701221 

8 . 9402420 

3 . 933342 

6 

2. 848972 

4 . 0000020 

5 .4598221 

1 . 0919621 

4 .976662” 

6 

3.823662 

4. 2000020 

6 .6686321 

1 . 3337321 

6 .289062” 

6 

5.090502 

4. 4000020 

8. 1450921 

1. 6290221 

7. 938552” 

6 

6.731032 

4. 6000020 

9 . 9484321 

1. 9896921 

1 .001012” 

5 

8.848602 

4 . 8000020 

1 . 2151022 

2.4302121 

1 .260972” 

5 

1.157402 

5 . 0000020 

1 .4841322 

2. 9682621 

1 . 586982“ 

5 

1.50726 2 

5 . 2000020 

1 .8127222 

3. 6254521 

1 .995522“ 

5 

1.955312 

5 . 4000020 

2. 2140622 

4. 4281321 

2 . 507182“ 

5 

2.527882 

5 . 6000020 

2 . 7042622 

5 .4085321 

3 .147582 

5 

3. 258142 

5 . 8000020 

3. 3030022 

6 . 6059921 

3 . 948652“ 

5 

4.187862 

6 . 0000020 

4. 0342922 

8 .0685821 

4.950142 

5 

5  .  3695  52 

S  2+ SAVE 


9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

8 

8 

8 

8 

8 

7 

7 

7 

7 

7 

7 

7 

7 

6 

6 

6 

6 

6 

6 

6 


6 
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‘ 

' 


. 


■ 
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(5.1) 

VYEXACTLUlV 

V  M+D  YE X ACT  X 
[1]  M+D pO 

C  2 ]  M+(*X),*X 

V 

RATIO  X 

RATIO  OF  ERRORS 

E ( Rl) / E ( B  2)  E (R1 )  / E(R2 )  E(B2) /E(R2) 


J 

V 

y 

V 

Y 

V 

1 . 059ff3 

7 .581 El 

1 .  820ff2 

3  .  7  0  IE  2 

1 . 718ff~l 

4. 881ff0 

5  .  1 1 92? 2 

8. 209ffl 

2 . 075ff2 

4  0  93  5ff2 

4. 053ff~l 

6. OllffO 

3 . 648ff2 

8 . 777E1 

2. 329E2 

6  .  6  7  OE  2 

6. 3  8  5ff  ~1 

7 . 599ff0 

2 . 975ff2 

9  .  282ffl 

2 . 5  7  9ff  2 

9  .  2  2  5ff  2 

8 .  670ff“l 

9  .  9  3  9E 0 

2 . 594ff2 

9 . 728ffl 

2  .  8  1 9ff  2 

1 . 327E3 

1 . Q87ff0 

1  .  365ffl 

2. 3  5 1E2 

1 . 0 1 2E2 

3 . 048ff2 

2 o  0  5  2E 3 

1  .  2  9  62?  0 

2  .  0  2  7ff  1 

2 . 1 8  6E 2 

1 . 046ff2 

3 . 263ff2 

3 o  6  85ff3 

1 . 493ff0 

3 . 522ffl 

2 . 066ff2 

1 .  0  7  8E 2 

3 . 463ff2 

1 . 063ff4 

1 . 676ff0 

9 . 885ffl 

1  .  976ff2 

1  .  102ff2 

3 . 649ff2 

1. 8  85ff4 

1 . 846ff0 

1 .  7 1  Iff  2 

1 . 907ff2 

1 . 125ff2 

3 . 821ff2 

5  o  67  5ff3 

2. 004ff0 

5 . 045ff 1 

1 .  8  52ff2 

1 .  1 4  5ff 2 

3 . 979E2 

3. 567ff3 

2. 149ff0 

3  .  1 1 6ff  1 

1 .  801 E 2 

1. 1 6  2E2 

4 . 1  2  5ff  2 

2  o  7  0  9E 3 

2 . 282ff0 

2.  3  3  Off 1 

1 . 771E2 

1 . 178ff2 

4. 259ff2 

2  .  2  4  42?  3 

2 . 405ff0 

1 . 905ffl 

1 . 7  4  OE 2 

1 . 1 9  2E2 

4 . 383ff2 

1  .  9  5  3E 3 

2. 518ff0 

1  .  6  39ff 1 

1 . 715ff2 

1 .  2  0  42?  2 

4.497ff2 

1 o  7  55ff3 

2 . 623ff0 

1 . 4  5  7ff  1 

1  .  6  9  2E 2 

1.  2 1  6ff  2 

4.60  3E 2 

1  „  6 1 OE 3 

2. 7  2  OE 0 

1 . 325ffl 

1 . 673ff2 

1  .  2  2  6ff  2 

4. 701272 

1 . 501ff3 

2. 809ff0 

1 . 225ffl 

1  .  657ff2 

1 . 235ff2 

4. 792E2 

1 . 41 5ff  3 

2 . 8  9  2E 0 

1 . 146ffl 

1 . 642ff2 

1 .  243ff2 

4 . 8  7  62?  2 

1.  3  46ff  3 

2 . 970ff0 

1 . 083ffl 

1 . 6  2  9ff  2 

1 . 251ff2 

4 . 95  5ff2 

1 . 290ff3 

3 . 042ff0 

1 .  0  3  Iff  1 

1 . 617ff2 

1 . 258ff2 

5 . 028ff2 

1 . 2  4  2E3 

3 . 109ff0 

9  .  8  8  Off  0 

1 . 6  0  7ff  2 

1 .  2  6  42?  2 

5 . 097ff2 

1. 202E3 

3  .  172ff0 

9  .  5  1 3ff  0 

1  .  597ff2 

1. 270ff2 

5  .  162ff2 

1 . 168ff3 

3  . 231ff0 

9 . 197ff0 

1 . 5  89ff2 

1. 275ff2 

5 . 222E2 

1 . 1 3  8^3 

3 . 287ff0 

8.  9 2 3ff 0 

1  .  581ff2 

1  .  280ff2 

5  .  279ff2 

1 .  1 1 2E 3 

3 . 339E0 

8 . 683ff0 

1  .  5  7  4#  2 

1 . 2  8  5ff  2 

5  .  333E2 

1 . 089ff3 

3  .  3  8  8E 0 

8.47  2ff  0 

1 .  5  6  7ff  2 

1  .  2  89ff  2 

5 . 384ff2 

1  .  0  6  8E3 

3 . 435ff0 

8 . 2  8  3ff  0 

1 . 561ff2 

1 . 293E2 

5 . 432ff2 

1 .  0  5  OS’  3 

3 . 479ff0 

8 . 1 14ff 0 

1  .  556ff2 

1 . 297E2 

5 . 478ff2 

1 . 033E3 

3 . 521ff0 

7.96  2ff  0 

1 . 551ff2 

1 .  3  0 Iff  2 

5 . 521ff2 

1. 018ff3 

3 . 560ff0 

7  .  8 2 5ff 0 

4 


SIGDIG  6 


' 

' 

. 

■ 

82 


(5.2) 

VY1PBIMELQ1V 

V  V+X  Y1PRIME  Y 

Cl]  V+Yl 2] , ( -0. 5x J[ 1] )+l . 5xJ[2] 

V 

X 

0  0.2  2  4 

Y 

1  1 


USE  1  PI 
X  RK 1  Y 


X 

Y 

7 

|2(Y) | 

|2(7)| 

2  .  0  0  0  0  0E~ 1 

1 . 2214020 

2 . 4  4  2  8  02“  1 

2 . 758182“ 

6 

5  .  516352 

4. 0  0  0  0  0E~ 1 

1. 4918220 

2.98364 2  1 

6 . 737682“ 

6 

1 . 347542 

6 . 0  0  0  0  0E~ 1 

1 . 8221120 

3.644212  1 

1 . 234412 

5 

2.46882 2 

8 . 000002“l 

2. 2255220 

4.451042  1 

2 . 010282“ 

5 

4. 020572 

1  .  0  0  0  0  0E 0 

2. 7182520 

5 . 436502“l 

3 . 069202“ 

5 

6 . 138412 

1 . 2  0  0  0  0E 0 

3 . 3200720 

6 „ 640142“l 

4.498482“ 

5 

8 . 996952 

1 . 4000020 

4. 0551420 

8.11027 2  1 

6 .410182“ 

5 

1  .  282042 

i . 6000020 

4. 9529420 

9 . 905892“l 

8 .947892“ 

5 

1  .  789582 

1  .  8  0  0  0  0E 0 

6 . 0495220 

1 . 2099020 

1 . 229512 

4 

2 .459022 

2 . 0000020 

7 . 3888920 

1 .4777820 

1 .668582“ 

4 

3 . 337172 

S+SA VE 

VY2P2XM2[[]]V 

V  V+X  Y2PRIME  $ 

]  7-<-(  -  0.  5xJ[l]  )  +  l.  5  x  J[ 2  ] 

V 

USE 2  K2 

X  RK2  Y 

X 

Y 

7 

12(2)1 

12(7) | 

2 . 000002“l 

1 . 2214020 

2 .442792“! 

6 . 064832“ 

6 

1 .608362 

4 . 0  0  0  0  OE  1 

1. 4918120 

2 . 9  83612“l 

1 . 527472“ 

5 

3 .974852 

6.00000 E  1 

1 . 8220920 

3.644162  1 

2. 877332“ 

5 

7  .  361162 

8 . 0  0  0  0  OE  1 

2 . 2254920 

4 . 450962“! 

4.806162 

5 

1.21082 2 

1 . 0000020 

2. 7182120 

5 . 436382“l 

7 . 510072 

5 

1  .  865852 

1 . 2000020 

3 . 3200020 

6 . 639962“l 

1 . 124422“ 

4 

2  .  758422 

1 . 4000020 

4. 0550420 

8 . 110  0  02“ 1 

1 . 633952 

4 

3 .962342 

1 .6000020 

4 . 9528020 

9.905512  1 

2 . 322362 

4 

5  .  5725  12 

1 . 8000020 

6 . 0493220 

1 . 2098520 

3 . 244752 

4 

7 .710692 

2 . 0000020 

7 . 3886120 

1 . 4777120 

4 . 471962“ 

4 

1 .  0  53272 

7 

6 

6 

6 

6 

6 

5 

5 

5 

5 

6 

6 

6 

5 

5 

5 

5 

5 

5 

4 


Sl+SA VE 


' 


. 
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(5.2) 

USE  2  R 2 
X  RK  2  Y 

X 

2 . 0  0  00  0E~ 1 
4. 000002~1 
6  .  0  0  0  0  02_1 
8.00000 E  1 
1 . 0000020 

1 . 2000020 
1 . 4000020 
1 . 6000020 
1 . 8000020 
2 . 0000020 

S  2+ SAVE 


Y 

1  .  2214020 
1 . 4  918120 

1  .  8220920 
2 . 2254920 

2  .  7182120 

3 . 3200120 
4 . 0550520 
4. 9528120 
6 . 0493420 
7 . 3886420 


V 

2. 4427  92~1 
2.98361 2  1 
3 . 644172_1 
4. 4  5  09  72^1 
5.436402  1 

6. 6  3  99  92~1 
8 . 110052_1 
9.905572  1 
1 . 2098620 
1 .4777220 


12(7)1 

6 . 0  6  4  8  42~6 
1 . 5  0  9  9  22” 5 
2. 815  042~5 
4.658742  5 
7 . 219172~5 

1 . 0  7  2  7  32~4 
1 . 5481 92_4 
2 .  18  6  7  72~4 
3 . 0  3  7  9  32~4 
4. 1651 32_4 


|2(7)| 

1 . 457282"6 
3  .  5  8  8  2 42” 6 
6 . 6227  22~6 
1 . 085972“5 
1 .  6  6  8  6 42” 5 

2 . 4  6  0  3  32- 5 
3 . 525472_5 
4.946822  5 
6 . 8  3  04  52~ 5 
9.311982  5 


SIGDIG  4 

6 

VYEXACTl  D]V 

V  M+D  YE X ACT  X 

[1]  M+D p0 

[2]  M+(.*X),*X 

V 

RATIO  X 

RATIO  OF  ERRORS 

2(21)72(22)  2(21)72(22)  2(22)72(22) 


Y 

7 

Y 

7 

Y 

7 

0.4548 

0  .  343 

0.4548 

0 . 3785 

1 

1  .  104 

0 . 4411 

0.339 

0 . 4462 

0.3755 

1.012 

1 .108 

0.429 

0.3354 

0.4385 

0.3728 

1.022 

1.112 

0.4183 

0 . 3321 

0.4315 

0 . 3702 

1.  032 

1 .115 

0 . 4087 

0  .  329 

0.4251 

0.3679 

1 .  04 

1.118 

0.4001 

0.3262 

0.4193 

0.3657 

1  .  048 

1.121 

0.3923 

0.3236 

0.414 

0 .3637 

1.055 

1.124 

0.3853 

0.3211 

0.4092 

0.3618 

1.062 

1.126 

0. 3789 

0.3189 

0 . 4047 

0  .  36 

1 .  068 

1.129 

0. 3731 

0.3168 

0. 4006 

0.3584 

1.  074 

1.131 

4 


SIGDIG  6 


' 

‘ 


• 

■ 
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(5-3) 

vyiPfliA/ffcmv 

V  V+X  Y1PEIME  Y 

[1]  7<-J[2],(llxy[l])-lOxJ[2] 

V 

X 

0  0.2  2  4 

Y 

1  1 


USE  1  R 1 
X  RK 1  Y 


X 

Y 

V 

lff(y) 1 

lff(y)  I 

2  .  0  0  0  0 Off- 1 

1 . 2  214  OE 0 

2 . 44280ff"l 

2 . 7581 8ff-6 

5 . 51635ff~7 

4 . 0  0  0  0  0E  1 

1. 4918  2E 0 

2.9  836  4ff  1 

6 . 73768ff"6 

1 . 347  54ff~6 

6. 0  0  0  0  0E~ 1 

1. 82211ff0 

3 .64421 E~1 

1 . 23441ff“5 

2 . 4  6  8  8  2ff ~6 

8.00000 E  1 

2.  2  2  5  5  2E 0 

4 . 45104£~1 

2. 0102  8ff~5 

4. 0  2  0  5 7ff ~*6 

1 . OOOOOffO 

2  .  718  2  5ff  0 

s^seso#  1 

3  .  0  6  9  2  Off  5 

6. 1 3841ff~6 

1 . 20000ff0 

3 . 32007ff0 

6 . 64014£’"1 

4 . 49  848ff~5 

8 . 99695ff”6 

1 . 4  0  0  0  0E 0 

4.  0  5  51 4#  0 

8.11027#  1 

6 . 41018ff~5 

1  .  2  8  2  04ff ~5 

1 . 60000ff0 

4.  9  5  2  9  42?  0 

9. 90589#“l 

8  .  9  4  7  8  9ff  5 

1 . 78958ff"5 

1 .  80000ff0 

6  .  0  4  9  5  2E 0 

1  .  2  0  9  9  Off  0 

1 . 2  29  5  Iff ~4 

2 . 45902ff~5 

2 . OOOOOffO 

7 . 38889ff0 

1 . 47778ff0 

1 . 66858ff  4 

3 . 337 1 7ff~5 

£*-£,4  VE 

V  Y 2PRIMEIU ] V 

V  y*-*  Y2PRIME  Y  ' 

1]  7«-(  llxj[  1]  ) -10xJ[  2] 

V 

USE2  K2 

X  RK  2  Y 

X 

Y 

V 

Iff(J) 1 

lff(y)  1 

0.2 

1 „ 22115 

0  .  244809 

0.000253851 

0. 000528847 

0 . 4 

1 . 49142 

0  .  299185 

0  .  000409276 

0 .000819951 

0.6 

1.82157 

0  .  36548 

0.000550561 

0 .00105619 

0  .  8 

2  .  22483 

0. 446412 

0. 000710759 

0. 00130404 

1 

2.  71737 

0 . 545249 

0 . 000907014 

0.00159221 

1.  2 

3. 31896 

0.665962 

0.00115269 

0. 00193815 

1 . 4 

4. 05374 

0.813397 

0 . 00146182 

0 .00235726 

1.  6 

4.95118 

0 . 993473 

0.00185102 

0  .  00286628 

1  .  8 

6 . 04731 

1 .21341 

0.00234081 

0 .00348493 

2 

7.3861 

1 . 48205 

0. 0029567 

0 . 00423694 

Sl+SA  VE 


. 
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(5.3) 


USE  2  R 2 
X  RK  2  Y 


X 

Y 

V 

|£( Y) | 

IE(V) | 

0.2 

1.22115 

0  .  244838 

0 .000253851 

0 . 000557316 

0.4 

1 . 49143 

0  .  299239 

0.000398924 

0. 000874106 

0 . 6 

1.8216 

0 .3655  58 

0.00051898 

0. 00113462 

0/8 

2 . 22489 

0. 446517 

0 . 000645967 

0 .00140888 

1 

2 . 71749 

0. 545385 

0 . 000794482 

0. 00172855 

1 .  2 

3 . 31914 

0. 666137 

0. 000973878 

0. 00211363 

1.4 

4. 05401 

0. 813622 

0 .00119265 

0. 00258203 

1 . 6 

4.95157 

0.99376 

0 . 00146017 

0 .00315337 

1 .  8 

6 .04786 

1 . 21378 

0.00178755 

0.00385085 

2 

7  .  38687 

1.48251 

0.00218827 

0 .00470249 

S2+-SA  VE 
SIGDIG  4 

6 

VYEXACTZUlV 

V  M+D  YEXACT  X 

[1]  M+D pO 

[2]  M+(*X),*X 

V 

RATIO  X 

RATIO  OF  ERRORS 

E(R1)/E(K2)  E(R1)/E{R2)  E(K2)/E(R2) 


Y 

V 

Y 

V 

Y 

V 

0. 01087 

0. 001043 

0. 01087 

0 .0009898 

1 

0 . 9489 

0 . 01646 

0  .  001643 

0.01689 

0 .001542 

1 . 026 

0.938 

0. 02242 

0  .  002337 

0 .02379 

0.002176 

1.061 

0.9309 

0.02828 

0.003083 

0.03112 

0 . 002854 

1 .  1 

0 .9256 

0  .  03384 

0.003855 

0.03863 

0 .003551 

1.142 

0. 9211 

0.  03903 

0 . 004642 

0. 04619 

0.004257 

1.184 

0.917 

0.  04385 

0 . 005439 

0. 05375 

0.004965 

1.226 

0.9129 

0  .  04834 

0 . 006244 

0.06128 

0 .005675 

1.  268 

0. 909 

0.05253 

0.007056 

0.06878 

0.006386 

1.31 

0.905 

0 . 05643 

0 . 007876 

0.07625 

0 .007097 

1.351 

0 . 901 

4 


SIGDIG  6 


. 

' 

. 

. 

86 


(5.4) 

VYIPRIMEIUIV 

V  V+X  Y1PRIME  Y 

[1]  7«-J[2],y[l] 

V 

X 

0  0.2  6  4 

Y 

1  1 

USE  1  R1 
X  RK1  Y 


X 

2 . 00000#~1 
4. 00000#_1 
6 . 0  0  0  0  0E~ 1 
8 . 00000#~1 
1  .  00000#0 

1 . 20000E0 
1 . 40000#0 
1 . 60000#0 

1  .  8  0  0  0  OE 0 

2. 0  0  0  0  OE 0 

2 . 20000E0 

2 .4000  OE 0 
2 . 60000#0 

2  .  8  0  0  0  OE 0 
3 .00000E0 

3  .  2  0  0  0  OE 0 

3 .4000  OE 0 
3 . 60000#0 

3 . 8  0  0  0  OE 0 
4.0000  OE 0 

4 . 2  0  0  0  OE 0 

4 . 4  0  0  0  OE 0 
4 . 6  0  0  0  OE 0 
4 . 80000#0 
5 . 0  0  0  0  OS’  0 

5  .  2  0  0  0  OE 0 

5 . 4  0  0  0  OE 0 
5 . 60000#0 
5 . 80000#0 
6 . 00000#0 


Y 

8 . 18733#“l 
6 . 70324#_1 
5 . 48817#_1 

4. 4  9  3  3  5#_1 
3.67885#  1 

3.  012  0  0E~ 1 
2 . 46602#“l 
2. 01902#_1 
1.  65304#_1 
1  .  3  5  3  4  0E~1 

1  .  10807#“l 
9.0721 4#  2 

7. 4  2  7  6  6#~2 

6.08 121 E  2 

4. 9  7  8  9 4#~2 

4. 0  7  6  4 3#~2 
3 .  3  375  1#_2 
2. 7  3  2  5  3#~2 

2  .  23721#_2 
i . 83168#  2 

1. 49966#”2 
1 . 22782#_2 
1 . 00526#_2 
8 . 23031 E  3 
6. 73848#~3 

5  .  517  0  2#~  3 
4 . 5169  IE _3 
3.69819#  3 
3 . 0  2  7  8  3#~  3 
2 . 47899#  3 


V 

1  .  6  3  7  4  7#~  1 

1  .  34065#~1 
1. 09763#~1 
8 . 98669#~2 
7  .  3  5  7  7  OE  2 

6 . 0  2  4  0  0E~  2 
4. 93205#~2 
4.038  03E  2 
3.  3  0  6  0  7#~2 
2 . 70679 E  2 

2  .  2161  4#~2 

1 .  8  144  3#  2 

1 . 4  8  5  5  3#~2 
1.21625 E  2 
9  .  9  5  7  8  8E~  3 

8 . 15285#~3 
6 .  675  01#_3 
5.46505#  3 
4. 47442#_3 
3.66336 E  3 

2  .  99931#~3 

2 .  45564#_3 
2 . 0 1 0  5 1#  3 
1  .  64607#~3 
1.34770#  3 

1  .  10  3  4  0#” 3 
9.03393#  4 
7  .  3  96  3  8#~4 
6  .  0  5  5  6  6#~4 
4.95797#  4 


l#m I 

2  .  5  8  0  2  6#~6 
4.22509#  6 
5  .  18  8  8  2#~6 
5 .  66434#~6 
5 .79697#  6 

5 . 69540#"6 
5 . 4401 7#~6 
5.09034#  6 
4. 68857#~6 
4. 2652 1#” 6 

3  .  8  412  7#-6 
3 . 43087#^6 
3,04305#  6 
2.  68  309#~6 
2.35364#  6 

2 . 0  5  547#~6 
1 . 7  8  8  0  6#_  6 
1  .  5  5  00  5#~6 
1 . 33958#  6 

1  .  1  5448#**6 

9 . 9  2  4  7  4#~  7 
8 . 51264#_7 
7.28637#  7 
6. 22496#~7 
5.30893#  7 

4 . 52045#~7 
3 . 84339#_7 
3 . 26325#_7 
2.76715#  7 

2  .  3  4  3  6  7#_7 


\E(V) | 

5 . 16052#"7 
8 . 45017#~7 
1.03776#  6 
1  .  13287#"6 
1 . 15939#”6 

1 . 13908#^6 
1.08803#  6 
1 . 01807#“6 
9. 37715#_7 
8.53041#  7 

7 . 68253#“7 
6.  8  617  5#~7 
6.08609#  7 

5  .  3  6  61 8#~7 

4. 7  0  7  2  8#~  7 

4. 110  9  4#~  7 
3 . 57611#"7 
3 . 10011#"7 
2 . 679 17#"7 
2 . 30897#"7 

1 .  98495#~7 

1 . 7  0  2  5  3#” 7 
1 . 45727#~7 
1  .  24499#“7 
1 . 06179#"7 

9. 0  4  09 1#~  8 
7 . 6  8  6  7  7#” 8 

6  .  5  2  6  5  0#~  8 
5 . 5  3429#~8 
4 . 68735#"8 


S+SA VE 


' 
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(5.4) 

V  Y 2PRIMEt □ ] V 

V  V+X  Y2 PRIME  Y 
[1]  V+Yl 1] 

V 

USE2  B 2 
X  RK2  Y 


X 

2.  0  0  0  0  0#_  1 
4 . 0  0  0  0  OS’”  1 
6 . 00000#_1 
8.00000#  i 
1 . 00000#0 

1  .  2  0  0  0  OE 0 
1 . 4000020 

1 . 6  0  0  0  OE 0 
1 . 80000#0 
2 . 00000#0 

2  .  2  0  0  0  OE 0 
2 . 40000#0 
2 . 6 0000# 0 

2  .  8  0  0  0  OE 0 
3 . 00000#0 

3  .  2  0  0  0  OE 0 

3 . 4  0  0  0  OE 0 
3 . 60000#0 
3  .  8  0  0  0  OE 0 
4 . 0  0  0  0  OE 0 

4 . 2  0  0  0  OE 0 
4 . 40000#0 
4 . 60000#0 
4.8000  OE 0 
5  .  0  0  0  0  OE 0 

5  .  2  0  0  0  OE 0 
5 .4000  OE 0 
5 . 60000#0 
5 . 80000#0 
6 . 00000#0 


Y 

8 . 1873 1E~1 
6  .  7  0  3  2  OE 

5  o  48812#_1 
4 . 49329#  1 
3 . 67880#~1 

3  .  0119  4#~  1 
2. 46597#_1 
2. 0189 1E~ 1 
1 . 65299#_1 
1. 35336#  1 

1  .  1  0  8  0  4#~  1 
9  .  07184#_2 

7 .  42742#~2 

6  o  08108#  2 

4.  9  7  8  8  0E~  2 

4. 0  7  6  3  3E~  2 
3  .  33746#_2 

2. 7  32 .5  4#  ”2 

2.  2  3  7  2  3E  J2 
1  .  83182#  2 

1. 49987#~2 
1.22811#  2 
1.  0  0  5  6 4#~2 

8.  23538#_3 

6. 7  4  4  8  2E  3 

5  .  5  2  4  9  6#~  3 
4. 52684#_3 

3.  71040#_3 
3.04286#  3 
2.49745#  3 


V 

1 . 63746#"1 
1 . 34064#"1 

1  .  0  9  7  6  2#~ 1 
8  .  9  8  6  5  QE~  2 

7  .  3  5  7  5  9#~  2 

6 . 02388#~2 
4. 9  319  4#” 2 
4.03793#  2 
3  .  3  0  5  9 1 E~  2 
2 . 70670#"2 

2  .  216  0 5#~2 
1. 81435#“2 
1 . 48546#_2 
1. 21619#  2 
9 . 95722#“3 

8  o  15221#"3 
6 . 67437#~3 

5 . 4  6  4  4  0#~  3 

4. 4  7  3  7  4#~3 
3.66262#  3 

2  .  9  9  8  5  0#” 3 
2.45471#  3 
2 . 00944#“3 
1 . 64482#~3 
1.34621#  3 

1. 1016  3#~  3 
9.01263#  4 
7. 37065#“4 
6 . 0  244  8#_4 
4. 9201 0#_4 


l#U)  I 

2.4775  7#~9 
1  .  1  5  7  3  5#~8 
2.64711#  8 
4  .  6  8  0  0  8#~  8 
7.25865#  8 

1 . 0421 4#- 7 
1 . 42423#“7 
1  .  8  8  3 1 1#~ 7 

2. 4  3  3  6  4#~  7 
3.09500#  7 

3 . 8  913  5#- 7 

4. 8  5  2  7  0#~7 
6 . 01603#_7 
7.42  67  2E~ 1 
9  .  1403  0#~7 

1  .  1 2  247#~6 

1  .  3  7  6  2  8#~  6 
1.68559#  6 

2  .  0  6  2  7  7#~6 
2.52291#  6 

3 .  0  8  4  4  6#~  6 
3 . 7699 1#~6 
4 . 60676#"6 
5 . 62857#_6 
6.87634#  6 

8 . 40014#~6 
1.02611#  5 
1  .  2  5  33  9#~5 
1  .  5  3  09  8#~5 
1  .  8  7  0  0  2E~  5 


I  #(7)  I 

6. 9711 5#-9 
1  .  2  918  3#~  8 
1. 83165#_8 
2.35758#  8 
2 . 9  0  6  5  7#~  8 

3  .  513  6  4#” 8 
4. 213  7  6#~  8 
5 . 04372#_8 
6 . 043  95#~8 
7 . 2  6  04  2#” 8 

8 . 7  46  7  3#~8 
1 . 05664#"7 
1  .  2  79  5  8#~7 

1  .  5  5  2  6  8#~  7 
1. 8871 3#~7 

2  .  29  6  5  2#~7 
2 .  79743#~7 
3.41008#  7 
4. 15  91 6#~  7 
5.07480#  7 

6 . 19382#"7 
7 . 56118#~7 
9  .  2  3 1 8  0#~7 

1 .  12728#~6 
1.37660#  6 

1 . 6811 6#~6 
2.05319#  6 

2.  5  0  7  6 1#~6 
3.0626  6#~6 

3 . 7  4  0  6  2#~6 


Sl+SA VE 


0 


■ 


* 
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' 
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(5.4) 

USE 2  R2 
X  RK  2  Y 


X 

2. 00000 2™1 
4 . 00000 2™1 
6  .  0  0  0  0  OF-! 
8 . 0  0  0  0  QE~ 1 
1 . 0000020 

1  .  2  0  0  0  OE 0 

1 . 4  0  0  0  OE 0 
1  .  6  0  0  0  OF 0 
1  .  8  0  0  0  OE 0 
2 . 0000020 

2 . 20000E0 

2 . 4  0  0  0  OF 0 
2.  6  0  0  0  OF  0 
2 .  8  0  0  0  OF 0 
3 . 0000020 

3 . 2000020 
3 ,4000  OF  0 
3  .  6  0  0  0  OF  0 
3  .  8  0  0  0  OF 0 
4. 0  0  0  0  OF 0 

4. 2000020 

4. 4000020 

4. 6000020 

4 . 8000020 
5 . 0000020 

5 .2000020 

5 .4000020 

5 . 6000020 

5 . 8000020 
6 . 0000020 


Y 

8 . 1873 12~1 
6 . 7  0  3  2  02^1 
5.488122  1 
4. 49  3  292~1 
3.678802  1 

3 . 011952~1 
2 . 465982_1 
2 . 018  9  72^ 1 
1 . 6  5  3  0  02” 1 
1 . 353372  1 

1 . 108052~1 
9 . 072012_2 
7 . 427622_2 
6 . 0  813  32^2 
4.979112  2 

4 . 076712~2 
3.337932  2 
2 . 7331 12~2 
2 . 237982  2 
1 . 832672~2 

1 . 5  0  0  9  02~2 

1 . 2  2  9  3  82^2 
1 . 0071 92^2 
8 . 2543 12_ 3 
6.767952  3 

5 . 5  5  3  2  22~  3 
4. 561352_3 
3.  7  5  2  5  52^3 
3 „ 094362~3 
2.560342  3 


V 

1. 6  3  7  4  62~1 
1. 34  0  6  42~1 
1.097622  1 
8. 986572~2 
7. 357582"2 

6. 023872"2 
4. 9  319  2F~  2 
4. 0379 12~2 
3. 305952_2 
2.706672  2 

2. 216  0  32” 2 
1 . 814312_2 
1.485422  2 
1. 2161 32~2 
9.956592  3 

8 . 151 442~3 

6.  6  7  3  4  32^3 
5 . 463262_3 
4. 4  7  2  3  42^3 
3.660912  3 

2. 9964 12” 3 
2 . 452172_3 
2.006342  3 
1. 641032™ 3 

1 . 3  415  82~  3 

1. 0  9  5  9  82~3 
8. 943582_4 

7.  2  8  6  3  22^4 
5 .  9  21 4  82_4 
4. 794302  4 


|2( Y) | 

2 . 026372"8 
6.  5  9  7  692” 8 
1.345092  7 

2 . 2  4  9  5  32” 7 
3.379462  7 

4. 7  5  5  7 12~  7 
6 . 413  4  62” 7 
8 . 402812~7 
1 . 079022_6 

1 . 3  6  6  032~6 

1 . 711 942~6 
2.  129  8  82~6 
2 . 636042™6 
3 . 250232~6 
3. 996692™6 

4. 9  0  5  0  52~6 
6 . Oil 492~6 
7 . 360192_6 
9 . 005122  6 
1 . 1 01 2 12~5 

1 . 346162"5 
1.645182  5 
2 . 010272~5 
2 . 456072™5 
3 . 000462™5 

3. 6  65  292~5 
4.477232  5 
5 . 468882™5 
6.6800  32~  5 
8 . 1 5  9272~5 


12(F) | 

2. 8  5  25  52~8 
5 . 326872"8 
7 . 6 1 1 5  52_ 8 

9 . 8  7 1442~8 
1 . 2257 12~7 

1 . 491322"7 
1 .  7  9  8  6  72~  7 
2. 16  3  3  62~7 
2.602722  7 
3  .  13  6  5  92~  7 

3 . 788192"7 
4. 5  8  5 1 72~  7 
5  .  5  6  0  7  02~ 7 

6. 7  549  32~7 
8.216602  7 

1 . 000502“6 
1 . 219262™6 
1 . 486742™6 
1  .  8137  32~6 
2. 213  3  72_6 

2. 7  017  42~  6 
3  .  2  9  8  4 42~6 
4. 027452_6 
4. 9 1 8  042_6 
6 . 005952  6 

7 . 334862“6 
8 . 9581 12~6 
1.094092  5 
1 . 336272_5 
1 . 632082  5 


S2+-SA  VE 


6 


SIGDIG  4 


■ 
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(5.4) 

1YEXACT\_U~]1 

V  M+D  YEXACT  X 

[1]  M+DpO 

[2]  M+(*-X) t-*-X 

V 

RATIO  X 

RATIO  OF  ERRORS 

E  (Rl)  /  E(B2)  E(R1) /E(R2)  E  ( B  2  ) /ff(  R2  ) 


I 

V 

Y 

V 

Y 

V 

1 .  041ff3 

7 . 403ffl 

1 . 273ff2 

i . 809ffi 

1 . 223 E~ 

1 

2.444 ff 

3  .  6 5  lff2 

6. 541ffl 

6 . 4  04ffl 

1  .  5  8  6ffl 

1.754 E 

1 

2 . 425ff 

1 . 960ff2 

5 . 666ffl 

3 . 858ffi 

i . 363ffi 

1 . 968ff~ 

1 

2 . 4  0  6ff 

1. 210E2 

4. 805 El 

2. 518ffl 

1 . 148ffl 

2 . 080ff~ 

1 

2.  3  8  8ff 

7 . 986ffl 

3 . 989ffl 

1 . 715ffl 

9 . 4  5  9E 0 

2 . 148ff" 

1 

2.37  Iff 

5 . 465ffl 

3 . 242ffl 

1. 198ffl 

7 .638ffo 

2. 191ff“ 

1 

2.35  6ff 

3 . 820ffl 

2  .  5  8  2E1 

8 . 482ff0 

6  .  0 4 92? 0 

2 . 221ff" 

1 

2.343 ff 

2. 703ffl 

2 . 018ffl 

6  .  0  5  8E 0 

4. 706ff0 

2.  2 4  Iff 

1 

2 . 3  3  Iff 

1 . 927ffl 

1. 551ffl 

4 . 345ff0 

3.  603ff0 

2 . 255ff“ 

1 

2. 3  2  2ff 

1 . 378ffl 

1 . 175ffl 

3.12  2E 0 

2 .  720ff0 

2.  2  6  6ff~ 

1 

2.31 5ff 

9 . 8  7 1E0 

8 . 783ff0 

2 . 244ff0 

2. 028ff0 

2  .  2  7  3ff- 

1 

2  .  3  09ff 

7 . 0  7  0E 0 

6 . 494ff0 

1 . 611ff0 

1 .497ff0 

2.27  8ff 

1 

2 . 304ff 

5 . 058ff0 

4. 7  5  6E 0 

1  .  154ff0 

1 . 0  9  42?  0 

2.282 ff~ 

1 

2. 3  0  Iff 

3 . 6i3ffo 

3 . 456ff0 

8.255 E  1 

7 .  944ff“l 

2.285 ff" 

1 

2.299 ff 

2  .  5  75ff0 

2. 4  9  42?  0 

5 . 889ff~l 

5 . 729ff~l 

2. 231 E 

1 

2. 297ff 

1 . 831ff0 

1  .  790ff0 

4. 1 9 l#- 1 

'  4 . 109ff“l 

2 . 233E~ 

1 

2 . 295ff 

1  .  299ff0 

1. 278ff0 

2 .  974ff"l 

2.933 E  1 

2  .  2  8  9ff  ~ 

1 

2 . 2  9  4ff 

9 . 19 6E~1 

9 . 0  9 1E~ 1 

2.  1  0  6ff  ~  1 

2 . 085ff“l 

2 . 2  9  Off  ~ 

1 

2.29  4ff 

6.494 ff“l 

6.442 E  1 

1.48  8E  1 

1. 477ff“l 

2.  2  9  Iff- 

1 

2  .  2  9  3ff 

4. 576 E  1 

4.550 E  1 

1. 048ff“l 

1 . 043ff“l 

2.29 lff- 

1 

2  .  29  3ff 

3 . 218ff“l 

3  .  2  0  5E~ 1 

7  .  3  7  3E~  2 

7 . 347ff-2 

2.  291ff" 

1 

2.293 ff 

2.258 E  1 

2 . 252ff~l 

5.17  4#  2 

5 .  162ff  2 

2. 2  9  Iff 

1 

2.292 ff 

1.5  8  2E  1 

1  .  5  7  9  ff  ~  1 

3.625 E  2 

3 . 6i8ff  2 

2.292 ff 

1 

2  .  29  2ff 

1 . 1 0  6ff  ~  1 

1 .  1  04ff  1 

2. 5  3  5ff  ~2 

2 .  531ff  2 

2. 292ff“ 

1 

2  .  2  9  2ff 

7 . 1 21E~  2 

7  . 713ff~2 

1  .  7  69ff  ~2 

1  . 768ff~2 

2 .292 ff" 

1 

2.292 ff 

5  .  3  8  Iff  ~2 

5  .  3  7  8E~  2 

1  .  23  3E~  2 

1  .  2  3  3E~ 2 

2 . 292ff" 

1 

2.292 ff 

3 . 746ff"2 

3.744 ff  2 

8 .  584ff_3 

8.581 E  3 

2.292 ff 

1 

2.29  2ff 

2.  6  0  4ff~2 

2  .  6  0  3E~2 

5  .  96  7E~  3 

5 . 965ff“3 

2  .  2  9  2ff 

1 

2.292 ff 

1 .  807ff~2 

1  .  8  0  7ff ~2 

4.142 ff  3 

4.14  2E  3 

2  .  2  9  2ff  ~ 

1 

2.292 ff 

1  .  2  5  3E~  2 

1 . 253E~2 

2.872 ff  3 

2.872 E  3 

2.29  2ff~ 

1 

2  .  2  9  2ff 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 


4 


SIGDIG  6 


'  . 


' 


\ 


. 


90 


(5-5) 

vyip#iAf#[[]]v 

V  V*-X  Y1PRIME  Y 

[1]  F^-J[2],(-0.5xJ[l])-i.5xy[2] 

V 

X 

0  0.  2  2  4 

Y 

1  “l 


USE  1  PI 
X  RK 1  Y 


X 

Y 

V 

i#m  i 

\E(V)  | 

2 . 00000#“l 

8. 18733#“ 

1 

“1. 63747#“ 

1 

2. 58026#“ 

6 

5.16052# 

4. 00000#  1 

6 . 70324# 

1 

1. 34065# 

1 

4. 22509#“ 

6 

8 .45017# 

6 . 00000#  1 

5 .48817#“ 

1 

1. 09763# 

1 

5 . 18882#“ 

6 

1 . 03776# 

8 . 00000#“l 

4.49335# 

1 

“8 . 98669#“ 

2 

5 . 66434# 

6 

1 . 13287# 

1 . 00000#0 

3 . 67885# 

1 

7. 35770# 

2 

5 . 79697#“ 

6 

1 .15939# 

1  .  2  0  0  0  0#0 

3 .01200#“ 

1 

“6 . 02400#“ 

2 

5 . 69540#“ 

6 

1 . 13908# 

1 . 4  0  0  0  0#0 

2. 46602#“ 

1 

“4.93205# 

2 

5 . 44017#“ 

6 

1 . 08803# 

1 . 6  0  0  0  OP  0 

2 . 01902#“ 

1 

4. 03803# 

2 

5 .09034#“ 

6 

1 .01807# 

1 . 8  0  0  0  0#0 

1 . 65304# 

1 

3.30607# 

2 

4. 68857#“ 

6 

9 . 37715# 

2 . 00000#0 

1 . 35340# 

1 

2. 70679# 

2 

4. 26521#“ 

6 

8 .53041# 

S+SA VE 

VY2PRIMEIU1V 

V  V+X  Y2PRIME  Y' 

]  7^-(-0.5xy[l])-l.5xy[2] 

V 

USE  2  K 2 

X  RK  2  Y 

X 

Y 

V 

l#U)  1 

l#( V) | 

2 . 00000#“l 

8. 18737#“ 

1 

“1 . 63748#“ 

1 

5 . 94026#“ 

6 

1 .60344# 

4. 00000#“l 

6 . 70329# 

1 

1 . 34067#“ 

1 

9  .  36932# 

6 

2.58982# 

6  .  0  0  0  0  OE  1 

5.48823# 

1 

“1. 09765#“ 

1 

1 .10364#“ 

5 

3  .  13356# 

8 . 00000#“l 

4. 49340# 

1 

8 . 98692#“ 

2 

1 . 14985#“ 

5 

3  .  36581# 

1 . 00000#0 

3 . 67891# 

1 

“7. 35793#“ 

2 

1.11654# 

5 

3  .  38440# 

1 . 2  0  0  0  OE 0 

3 . 01205#“ 

1 

“6 . 02421#“ 

2 

1 . 03355#“ 

5 

3 .26168# 

1 . 40000#0 

2 . 46606#“ 

1 

“4. 93224#“ 

2 

9 .22270#“ 

6 

3  .  05056# 

1 . 60000#0 

2 . 01904#“ 

1 

4.03821# 

2 

7  .  97741# 

6 

2.78918# 

1  .  8  0  0  0  OE 0 

1 . 65306# 

1 

“3. 30623#" 

2 

6 . 70293#“ 

6 

2 . 50456# 

2 . 00000#0 

1 . 35341#“ 

1 

“2. 70693# 

2 

5 . 46765#“ 

6 

2. 21541# 

7 

7 

6 

6 

6 

6 

6 

6 

7 

7 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 


Sl+SA  VE 


‘ 


t  >50  '  i  . 


■ 


£  . 


.3. 


. 


' 


91 


(9.5) 

USE 2  /?2 
X  RK  2  Y 


Z 

2 . 00000£"l 
4 . OOOOOF"! 
6  .  0  0  0  0  OS'”  1 
8 . 0  0  0  0  QE  1 
1 . OOOOOtfO 

1 . 20000E0 
1 .4000  OE 0 
1 . 60000E0 
1 . 80000E0 
2. OOOOOffO 


y 

8  .  18  7  3  7E~1 
6 . 70330 E~1 
5 . 48823F“l 
4 . 49  34  IE  1 
3 . 67891E  1 

3  .  012 0 62?”  1 
2. 4  6  6  0  7E~ 1 
2. oigoeE"! 
1  .  65307J?  1 
1  .  3  5  3  4  2E~ 1 


V 

1 .  63748£’“l 
1.  34066£’~1 
1 . 09765£~1 
8 . 98689£"2 
7  .  35790£’"2 

6 . 02419F"2 
4. 9  3222E~  2 
4. 03819£’"2 
3. 30621£“2 
2.70692 E  2 


iff(y) l 

5 . 94026£’~6 
9  .  50904£’"6 
1 .  13  91 SE~ 5 
1 . 210  0  9E~ 5 
1 . 2  017  2E~  5 

1 .  1 4  2  0  251”  5 
1  .  0  512  5E~5 
9 . 4  3  9  0 9#” 6 
8  .  30116£’_6 
7.16769E’  6 


I  E(V)  | 

1 . 44117£~6 
2 . 33808£"6 
2.84276 E  6 
3 . 06981£"6 
3  .  10  5  0  0E~6 

3 . 01195£’"6 
2.  8  3  7 4 0£’"*6 
2.  615  2  0E~6 
2. 36950£"6 
2  .  1171  6£’-6 


S  2+ SAVE 


SIGDIG  4 

6 

VYEXACTLUlV 

V  M+D  YEXACT  X 

[1]  M+D pO 

[2]  M+{ *-X) , -*-X 

V 


RATIO  X 

RATIO  OF  ERRORS 

E(R1) /E(K2)  E(R1) /E(R2)  E(K2 ) /E(R2 ) 


y 

V 

y 

V 

y 

V 

0.4344 

0.3218 

0 . 4344 

0.3581 

1 

1  .113 

0.4509 

0. 3263 

0 . 4443 

0. 3614 

0.9853 

1.108 

0.4702 

0. 3312 

0.4555 

0.3651 

0.9688 

1 .102 

0.4926 

0. 3366 

0.4681 

0  .  369 

0.9502 

1.096 

0.5192 

0 . 3426 

0.4824 

0 . 3734 

0.9291 

1 . 09 

0  .  5511 

0. 3492 

0.4987 

0.3782 

0.905 

1 .083 

0 . 5899 

0. 3567 

0.5175 

0.3835 

0 . 8773 

1 .075 

0 .6381 

0.  365 

0.5393 

0.3893 

0 . 8451 

1 . 067 

0 . 699  5 

0  .  3744 

0 . 5648 

0.3957 

0. 8075 

1 . 057 

0 .7801 

0.385 

0.5951 

0 .4029 

0.7628 

1  .  046 

4 


SIGDIG  6 


- 


?  : 5T  i .  3 


4 

* 


- 


92 


(5.6) 

VYiPflimDlv 

V  V+X  Y1PRIME  Y 

[1]  V+Y12~]  ,(llxy[i]  )  +  10xY[2] 

V 

X 

0  0.2  2  4 

7 

1  "1 


USE  1  PI 
X  RK1  Y 


X 

Y 

V 

U?(7) | 

ItfOOl 

2 . 0  0  0  0  0E~ 1 

8  .  18733P"l 

“1. 63747£“1 

2.  5 8 0 2 6#"*6 

5  .  16052 E~ 

4  .  00000£_1 

6. 70324P_1 

“1. 34065£"l 

4. 22509E~6 

8 . 45017P" 

6 . 00000P“1 

5 . 48817P~1 

”1. 0  9  7  6 327~  1 

5  .  18882£’"*6 

1 .03776 E~ 

8.  0  0  0  0 OP” 1 

4.49335P  1 

” 8 . 9  8  6  6 9£’"~2 

5 . 66434P"6 

1. 13287 E~ 

1 . 00000F0 

3. 67885P_1 

~7 . 3  5  7  7  0E~  2 

5 . 7  9  6  9  7E~  6 

1.15939 E~ 

1  .  2  0  0  0 0P0 

3. 012  0  OP-! 

"6 . 0  24  0 OP- 2 

5 . 69  5  40£’"*6 

1. 13908P" 

1 . 40000F0 

2.46602P  1 

"4. 93205P_2 

5. 44017£’”6 

1 . 08804P" 

1 .  60000P0 

2. 01902P  1 

~4 . 0  3  8  0  3E~  2 

5  .  09  033E~  6 

1. 01807P" 

1 .  80000P0 

1. 65304P“l 

”3.  30607£'~2 

4 . 6  8  8  5  5P~6 

9.37757 E~ 

2 .  00000P0 

1.35340P  1 

"2  .  7  0679£’~2 

4. 26504£~6 

8. 53399P" 

S+SA VE 

7Y2PRIMELU17 , 

V  V+X  Y2PRIME  Y 

]  y«-(llx7[l])  +  10xy[2] 

V 

USE 2  K2 

X  RK2  Y 

X 

Y 

V 

|P(7) 1 

\E(V) I 

2. 0  0  0  0 OP- 1 

8 . 19  01  OP- 1 

"1.  63140£’~1 

2  .  7 9  48  72?~4 

6. 06065P" 

4 . 0  0  0  0  OE  1 

6 .72882 E  1 

1. 28432P”! 

2. 5  6 1 6  IP- 3 

5 . 63219P" 

6 . 0  0  0  0  0E~1 

5 . 7  0  6  2  6P” 1 

"6. 16  7  7  2E~  2 

2.  1814  2E~~  2 

4. 80851P" 

8 . 00000P_1 

6  .  34042£’”l 

3.  17450ff“l 

1 . 84713P"! 

4. 07316P" 

1.  00000P0 

1 . 93125P0 

3. 37404P0 

1  .  56337P0 

3 . 44762P0 

1 . 2  0  0  0  OE 0 

1 . 35329P1 

2 . 91190P1 

1 . 3231 7 El 

2. 91792P1 

1 . 40000P0 

1 . 12233E2 

2. 46910P2 

1 . 11987P2 

2. 46959P2 

1 .  60000PO 

9 . 48006P2 

2. 09011P3 

9 . 47804P2 

2. 09015P3 

1  .  8  0  0  0  0P0 

8 . 02193P3 

1 . 76900P4 

8 . 02177P3 

1 . 76900F4 

2  .  00000P0 

6 . 78926P4 

1 . 49720P5 

6 . 78925P4 

1 . 49720P5 

7 

7 

6 

6 

6 

6 

6 

6 

7 

7 

4 

3 

2 

1 


Sl+SA  VE 


ivosref  .£ 


' 


.  ’ 


' 


:  1 


.  ■  ■  .  >  '  , 


‘ 


93 


(5.6) 

USE  2  R 2 
X  RK  2  Y 


X 

Y 

V 

\E(Y) | 

1^(7) 1 

2  .  0  0  0  0  02~  1 

8 . 1901 0E~ 1 

”1. 6  3  1  3  IS1”  1 

2. 79  4  8  7£’~4 

6 . 15040£~4 

4. 000002"l 

6 . 7291 0E~ 1 

”1. 2  8  3  6 42~1 

2.  59014£'~3 

5 . 69966£"3 

6 . 000002"l 

5 . 7  0  8  8  22"l 

”6  .  119  6  92~2 

2. 20701£”2 

4. 8  5  6  5  5E~  2 

8 . 0  0  0  0  0E~1 

6.  3  5  9  4 IE- 1 

3.  207742"l 

1. 86612ff“l 

4. 106  4  0E~1 

1. 0000020 

1 . 9445820 

3. 3959620 

1 . 5  767  0E0 

3 . 46954£,0 

1 . 2  0  0  0  OE 0 

1.  3622021 

2 . 9252221 

1 . 33208^1 

2 . 9  3124£’l 

1 . 4  0  0  0  OE 0 

1 . 1278622 

2. 47595F2 

1 . 12540F2 

2. 47644F2 

1 . 6000020 

9  .  5  09  8  722 

2 . 09217F3 

9 . 50785E2 

2. 09221£3 

1  .  8  0  0  0  0E0 

8. 0328123 

1. 76758£4 

8. 03264E3 

1 . 76759£,4 

2 .00000E0 

6  .  7  8  6  3  424 

1 . 49  3  3  3E5 

6. 78633F4 

1 . 4933425 

S  2+ SAVE 


SIGDIG  4 


VYEXACTZU1V 

V  M+D  YEXACT  X 

[1]  M+D pO 

[2]  M+(*-X) ,-*-X 

V 

RATIO  X 

RATIO  OF  ERRORS 

E(R1)/E(K2)  2(21)/2(i?2)  E(K2)/E(R2) 


X 

V 

J 

V 

Y 

V 

9 . 2322" 

3 

8. 5152" 

4 

9 . 2322" 

3 

8. 3912" 

4 

1 .00020 

9 . 8542"l 

1 . 6492" 

3 

1 . 5002" 

4 

1. 6312 

3 

1.4832" 

4 

9 . 8902"l 

9 .8822"! 

2. 3792" 

4 

2 .1582" 

5 

2.3512" 

4 

2. 1372" 

5 

9 . 8842"l 

9 . 9012"l 

3. 0672" 

5 

2. 7812" 

6 

3. 0352" 

5 

2.7592” 

6 

9. 8  9  82"l 

9 . 9192"1 

3 . 7082" 

6 

3. 3632" 

7 

3 .6772 

6 

3.  3422 

7 

9 . 9152”l 

9 . 9372"l 

4. 3042" 

7 

3 .9042" 

8 

4.2762” 

7 

3 .8862" 

8 

9. 9  3  32"l 

9 . 9552"l 

4. 8582" 

8 

4. 4062 

9 

4 . 8342” 

8 

4. 3942" 

9 

9.95 12" 1 

9 . 9722"l 

5 .3712" 

9 

4 .8712" 

10 

5. 3542" 

9 

4.8662” 

10 

9. 9  6  92" 1 

9 . 9902”1 

5 . 8452" 

10 

5 . 3012" 

11 

5 .837 2" 

10 

5 . 3052" 

11 

9 . 9862"1 

1.00120 

6 . 2822" 

11 

5 .7002" 

12 

6 .2852" 

11 

5. 7152" 

12 

1 .00020 

1 . 00320 

4 


SIGDIG  6 


■  .  ' 


. 

V'. 

■ 

■ 

‘ 

• 

. 


- 

. 


V 


,  * '  .  •' 


.  . 


■ 


r - 8  .  ^ 
•  V  .  c 
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(5.7) 

V  Y 1PRIMELU ] V 

V  V+X  Y1PRIME  Y 
[1]  V*-YZ  2  ]  ,  Y  [  1  ] 

V 

X 

0  0.2  2  4 

Y 

2  0 


USE  1  31 
X  RK 1  Y 


X 

Y 

7 

1  3  (  Y )  | 

1300  | 

2  .  0  0  0  0  0E~1 

2 . 0401330 

8 . 0  5  3  3  33~2 

1 .  779143" 

7 

1 .  067693" 

6 

4. 0  0  0  0 0E~ 1 

2. 1621430 

1.  64  2  9  93**1 

2 . 512593" 

6 

2  .  192553" 

6 

6.00000 E  1 

2 . 3709230 

2. 5  4  6  5  83"l 

7. 155303" 

6 

3  .  506593" 

6 

8 . 0  0  0  0  0E  1 

2. 6748630 

3.  5  5  23  73~1 

1 . 443853" 

5 

5 .153433 

6 

1. 0000030 

3 . 0861430 

4. 7  007  33~1 

2 .489513" 

5 

7.297803" 

6 

1 .  2  0  0  0  OS'  0 

3 . 6212730 

6.0377  43~1 

3  .  928943" 

5 

1 . 013603" 

5 

1 . 4  0  0  0  OE 0 

4. 3017430 

7.  617  0  73**1 

5 .866173" 

5 

1 .  390843" 

5 

1 . 6  0  0  0  OE 0 

5  .  1548430 

9 . 502083"! 

8 .438863" 

5 

1 . 891393" 

5 

1 . 8000030 

6 . 2148330 

1. 1768430 

1 .182623“ 

4 

2. 552793" 

5 

2.  0  0  0  0  0E0 

7 .5242330 

1. 4507130 

1 . 625933 

4 

3 .422473" 

5 

S+SA VE 

VY2PRIMELU1V 

V  V*-X  Y2PRIME  Y 

]  V+Y [1] 

V 

r 

USE  2  3  2 

Y  332  y 

y 

Y 

V 

1  3  (  Y )  | 

|3(7)| 

2  .  0  0  0  0  0E~ 1 

2. 0401330 

8 . 0  5  3  4  43~2 

1 .270403" 

10 

3 . 053273" 

10 

4 .  000003"l 

2. 1621430 

1 . 643013"l 

1  .  589773" 

9 

3 . 496583" 

9 

6 .  000003"l 

2. 3709330 

2.546613  1 

7 . 366483 

9 

9 . 813093" 

9 

8  .  0  0  0  0  03~ 1 

2. 6748730 

3.55242 3  1 

2. 077913" 

8 

1. 973943" 

8 

1 .  0000030 

3. 0861630 

4 . 700803"1 

4. 575083" 

8 

3.403403" 

8 

1 .2000030 

3. 6213130 

6 . 037843"l 

8 .709883" 

8 

5 . 377423“ 

8 

1 . 4000030 

4. 3018030 

7 . 617213"l 

1 . 508813" 

7 

8 .041953" 

8 

1 . 6000030 

5 . 1549330 

9 . 5  0  2  2  73~1 

2 . 448173" 

7 

1 . 158973" 

7 

1  .  8000030 

6 . 2149530 

1 .1768730 

3 .788023" 

7 

1 . 627143" 

7 

2. 0000030 

7. 5243930 

1 .4507430 

5 . 655523" 

7 

2 . 241013" 

7 

Sl+SA  VE 


■ 


'  a  ;  a  .  . . ' 


•  o  ■  . 


. 

■ 

95 


(5.7) 

USE 2  i?2 

y  i?X2  y 

y 

2. 0  0  0  0  0Z?~  1 
4 . 0  0  0  0  OS’  1 
6 . 0  0  0  0  0E~ 1 
8 . 00000£'"1 
1 . 00000£0 

1  .  2  0  0  0  OE 0 
1 . 4000020 
1  .  6  0  0  0  OF  0 
1  .  8  0  0  0  OE  0 
2 .OOOOOEO 

S2+SA VE 


Y 

2.  04013E0 
2. 1621 42? 0 

2. 3  7  0  9  3E 0 
2. 67487F0 
3. 086i6£0 

3  .  6  21  3 12?0 

4 . 3  018  OE 0 

5  .  15  49  32? 0 

6  .  2 1  4  9  5Z?0 
7 . 52439£0 


V 

8 . 05344£’"2 
1.643  01E  1 
2 . 54662F~1 
3.55242 E  1 
4. 70081£'"l 

6 . 03785£’“1 
7 . 617  21E~1 
9  .  5  0  2  2  7Z?- 1 
1  .  17687£’0 
1 . 4507420 


|F(y) | 

3  .  5  4  21 92?”  8 
9  .  8  4  5  3  02?~8 
1  .  875012”7 
3.029022  7 
4 . 4  6  8  0  72” 7 

6 . 231662“7 
8  .  3  7  81 52”  7 
i . 098662”6 
1 . 415  9  62” 6 
1 . 8  0  2  7  52”  6 


\E(V)  I 

2 . 703492“8 
5.0537  82  8 
7 . 241392”8 
9 . 435622~8 
1 . 179462”7 

1 . 447472“7 
1.763872  7 
2.  146  5  3E~  7 

2  .  615  7  6E~  7 

3  .  1  95  402”7 


SIGDIG  4 

6 

VYEXACT CD]V 

V  M+D  YEXACT  X 

[ 1 ]  M+D p  0 

[2]  AM  (*X)+*-X)  ,  (*X)-*-X 

V 


RATIO  X 


RATIO  OF  ERRORS 


E(R1) /E(B2) 
Y  V 


2(21)/2(22 ) 
I  V 


E(B2) /E(R2 ) 
y  V 


1 .40023 
1  .  58023 
9.71322 
6 . 94922 
5 . 44122 

4. 51122 
3  .  8  8  822 
3 . 44722 
3 . 122E2 
2 . 87522 


3. 49723 
6. 271E2 
3 . 57322 
2 . 61122 
2 . 14422 

1 . 88522 
1 . 72922 
1 .63222 
1  .  5 6 922 
1 . 52722 


5 . 02320 
2. 5  5  2E1 
3 . 81621 
4. 76721 
5  .  5  7  221 

6. 30521 
7. 00221 
7 .68121 
8 . 352 El 
9 . 01921 


3.  94921 
4 . 338 El 
4. 8  4  221 
5 .46221 
6. 18721 

7 .00321 
7 .88521 
8 .81121 
9 . 75921 
1 .07122 


3  .  5862_3 
1 . 6 1 52”  2 
3 . 9292~2 
6.860 E  2 
1 . 0242”l 

1 . 3982“l 
1 . 8012~1 
2. 228 E~1 
2 .675 2”l 
3.1 372  1 


1 .129 E  2 
6. 9 1 92”2 
1 . 3552"l 
2.09  22” 1 
2. 8862“l 

3 . 7152”l 
4 . 5592~i 
5 . 3992”l 
6. 2212“l 
7 . 0132”l 


4 


SIGDIG  6 


. 


\ '  . 

■ 


. 


96 


(5.8) 

VJ1P/?JM£[[]]V 

V  V+X  YIP RIME  Y 

[1]  V+Yl  2]  ,  7 C  1  ]  -  2  y-SIN  X 

V 


y 

0  0.2  2  4 

y 

1  0 


USE  1  R 1 
X  RK 1  Y 


X 

2 . 00000P~1 
4 . 00000P  1 
6. 00000P~1 
8 . 0  0  0  0  OE  1 
1  .  00000P0 

1  .  2  0  0  0  OP  0 
1 . 40000F0 
1. 60000P0 
1  .  8  0  0  0  OE 0 
2 . 00000P0 

S+SA VE 


Y 

1 . 01740P0 
1 . 05975P0 
1 . 113  4  6E 0 
1  .  16  6  7  OE 0 
1 . 20937P0 

1 . 23325F0 
1 . 23207P0 
1 . 20150P0 
1  .  13  91 827  0 
1 .  04467P0 


V 

3. 22667P~2 
5. 01475P_2 
5 . 53043P_2 
4. 94754P“2 
3. 44851P~2 

1 . 22341P"2 
1  .  5  3  2  3  5P~2 
4. 62155P  2 
7  .  849  5  1E~  2 
1.  10290P"! 


|p(y) I 

4.  27987£’“6 
7 . 62708£’~6 
1 . 03957P"5 
1  .  2  8  8  37E~  5 
1  .  5  3  5  0 OP- 5 

1 . 8  0  2  8  6P- 5 
2 . 1 1 4  05P~5 
2 . 49053P"5 
2 . 95515£~5 
3  .  5  3  2  7  2E~5 


[1] 


VY2PRIMEIU1V 
V  V+-X  Y  2  PR  I  ME  Y 
V+Yll]-2xSIN  X 


USE  2  B 2 
X  RK  2  Y 

X 

2 . OOOOOE'l 
4. 00000£"l 
6. 00000P"! 

8  .  00000P“l 
1 . 00000P0 

1 . 2  0  0  0  OE 0 
1 . 40000P0 
1.  6  0  0  0  OE 0 
1 .  8  0  0  0  OE 0 
2. 00000P0 


y 

1 . 01740P0 
1 . 05974P0 
1 . 11345P0 
1 . 16668P0 
1 . 20935P0 

1 . 23323P0 
1 . 23204P0 
1 . 20147P0 
1 . 13914P0 
1 . 04463P0 


V 

3.  2  2  6  7  2E~  2 
5 . 014  8  2E~  2 
5  .  5  3  0  4  QE~2 
4. 94755P"2 
3. 44844P"2 

1  .  2  2  3  25E~2 
1. 53263£~2 
4.62197P  2 
7  .  8  5  0  0  8E~  2 
1. 10297 E~1 


I E(Y) | 

3  .  89826£~7 
7 . 66000£“7 
1  .  13  3  9  2E~6 
1 . 49  9  7  8E~6 

1 .  87124£’"6 

2. 25807P~6 

2 .  67281P  6 
3 . 13143P~6 
3  .  6  5  4  01£’-6 
4. 26556£~6 


\E(V)  \ 

5 . 1018  6E~ 7 
6. 6  8  46  IP” 7 
5.26346P  7 
1  .  22739£~7 
5.1381 8E  1 

1 . 36379P"6 
2 . 41599£“6 
3.66757P  6 
5 . 12461£“6 
6. 8  0  3  2  2#” 6 


\E(V) | 

8.1717  7E~9 
4 . 3  6  5  0  QE~ 9 
3  .  713  9  2E~  8 
8  .  9  7  7  3  3£’~  8 

1  .  62039£’"7 

2  .  5  3  8  9  4£’~ 7 
3 . 65620P'7 
4 . 97945E'~7 
6.52219P  7 
8  .  3  0  6  2  2£’~7 


Sl+SA  VE 


' 


■ 


97 


(3.8) 

USE 2  R2 
X  RK  2  Y 


y 

2 . 00000£~1 
4.0000  0E  1 
6 . 00000E"1 
8 .00000E~1 
1 . 00000£0 

1 . 20000E0 
1 . 4  0  0  0  OE 0 
1 .6000020 
1  .  8  0  0  0  OE 0 
2 . 0  0  0  0  OE 0 

S2+SA VE 


Y 

1 .  0174020 
1. 0597420 
1. 1134520 
i. 1666820 
1. 2093520 

1.  2332320 
1. 2320420 
1. 2  014  62 0 
1.  13  91 420 
1.  0446220 


V 

3 .  226722~2 
5 . 014822~2 
5.530472  2 
4.9475  42~2 
3 . 448432~2 

1 .  223222~2 
1. 532672~2 
4. 622022~2 
7 . 8  5  01 62~ 2 
1 .  1  02982*"l 


12m  i 

9 . 411062~7 
1  .  8  343  82~6 
2. 690952~6 
3  .  523282~6 
4. 3467  92_6 

5 . 181472~6 
6 . 053362~6 
6  .  9  9  6 1 12~6 
8  .  052452~6 
9 . 275852_6 


I  E(V)  | 

3  .  3  2  0  0  8E~  8 
1 .  711282_8 
4 . 6  3  01 0E~  8 
1  .  553822_7 
3 .  088382"7 

5 . 058302_7 
7 . 462172~7 
1 .  030822_6 
1 . 361792_6 
1 . 743022_6 


SIGDIG  4 

6 

vy£,y^c!Z,[D]v 

V  M+D  YEXACT  X 

[1]  M+D  pO 

[2]  M*-(  (*- X)  +SIN  X)  A-*~X)+COS  X 

V 


RATIO  X 


RATIO  OF  ERRORS 


2(21)/2(S 2)  2(21)/2(22)  E(B2)/E(R2) 


y 

V 

y 

V 

y 

V 

10.98 

62 . 43 

4. 548 

15 . 37 

0.4142 

0. 2461 

9.957 

153.1 

4.158 

39 . 06 

0.4176 

0.2551 

9  .  168 

14. 17 

3.863 

11.37 

0.4214 

0.8021 

8.  59 

1.  367 

3.657 

0.7899 

0.4257 

0. 5778 

8 . 203 

3.171 

3.531 

1 .664 

0.4305 

0.5247 

7.984 

5 .371 

3 .479 

2.  696 

0 .4358 

0 .5019 

7  .  909 

6.  608 

3.492 

3.238 

0.4415 

0 . 49 

7.953 

7  .  365 

3.56 

3.558 

0 . 4476 

0.4831 

8 . 087 

7.857 

3.  67 

3.763 

0.4538 

0.4789 

8.  282 

8.  191 

3.809 

3.  903 

0.4599 

0.4765 

4 


SIGDIG  6 


' 


. 

■ 


98 


(5.9) 

VYl PR I ME lUll 

V  V+X  Y1PRIME  Y 
Cl]  7+Y[2],-Y[l] 

V 

X 

0  0.2  2  4 

Y 

0  1 


USE1  R 1 
X  RK1  Y 


X 

Y 

V 

1 E( Y) | 

•  152001 

2 . 000005"l 

1  .  9  8  6  6  75" 

1 

1 .96013 E~ 

1 

2  .  6641251" 

6 

1.776535 

4.00000 5"l 

3. 8941 3E~ 

1 

1 . 84212 E~ 

1 

5  .  18  6  7  352" 

6 

2 . 465325 

6.00000 E  1 

5  .  64  6  3  55 

1 

1 .65068 E~ 

1 

7  .  257675’" 

6 

6.715585 

8 . 0  0  0  0  0E~ 1 

7 . 173475 

1 

1  .  39343£’" 

1 

8  .  5945  652" 

6 

1 .262065 

1 . 0  0  0  0  0E 0 

8 . 414625" 

1 

1 . 08062£’" 

1 

8  .  9619852" 

6 

1 .973005 

1 . 2000020 

9 . 32031 5" 

1 

7 . 24743 E~ 

2 

8.  1881851" 

6 

2.74771 5 

1 . 4000050 

9 . 854445" 

1 

3 . 39969F" 

2 

6  .  17815 E~ 

6 

3 .521335 

1 . 6  0  0  0  0E 0 

9 . 9957 IE 

1 

"5 . 

3 

2  .  9  2  2  4 252" 

6 

4. 224715 

1  .  8  0  0  0  0E 0 

9  .  7  3  8  495" 

1 

4. 54356 E~ 

2 

1 . 4  9  8  9  652" 

6 

4.788695 

2 . 0000050 

9  .  093045" 

1 

“8.  3224252" 

2 

6  .  917  6 452" 

6 

5 .148535 

S+-SA  VE 

VY2PRIMEIU1V 

V  V+-X  Y2PRIME  Y 

]  V+-YI 1] 

V 

USE  2  B2 

X  RK2  Y 

X 

Y 

V 

|£(Y) 1 

|5(7)| 

2  .  000005"l 

1 . 986695" 

1 

1.  9  6  0 1 352" 

1 

2. 538275’" 

9 

7.098425 

4. 0  0  0  0  0E~1 

3  .  894185" 

1 

1 .  8  4  2 1 2£2~ 

1 

1 .201405’" 

8 

1.378285 

6.  0  0  0  0 05"  1 

5 . 64643 E 

1 

1  .  6  5  0 6  752~ 

1 

2. 7  7  8  0  552" 

8 

1 .950955 

8.000005  i 

7  .  17  3  5  65" 

1 

1 .  3  9  3 4 Iff" 

1 

4. 8  6  6  3 852" 

8 

2.379295 

1 . 0000050 

8 . 414  7  IE 

1 

1 . 0  8  0  6  052 

1 

7 .3048052 

8 

2.621845 

1  .  2  0  0  0  0E0 

9 . 320395" 

1 

7 . 24716 E~ 

2 

9  .  895  3652" 

8 

2.650045 

1 . 4000050 

9.8545  OE 

1 

3.  3993552 

2 

1 .241665’" 

7 

2.447075 

1 . 6  0  0  0  OE 0 

9 . 9  9  5  7 45" 

1 

"5.83988 E~ 

3 

1 . 463685’" 

7 

2.010605 

1 . 8000050 

9  .  7  3  8  4  8E~ 

1 

"4.  5440452" 

2 

1. 6  3  2  8  55’" 

7 

1 .353015 

2. 0000050 

9.09298 E 

1 

"8.  3  2  29  451 

2 

1 .7283152" 

7 

5 .011595 

8 

7 

7 

6 

6 

6 

6 

6 

6 

6 

9 

8 

8 

8 

8 

8 

8 

8 

8 

9 


Sl+SA  VE 


. 


* 


' 


99 


(5.9) 

USE 2  i?2 
X  RK2  Y 


X 

2 . 0  0  0  0  0E~ 1 
4.  0  0  0  0  OS’"! 

6  .  00000£_1 
8 . 0  0  0  0  OE  1 
1 . 0  0  0  0  OE 0 

1 . 2  0  0  0  OE 0 
1. 4  0  0  0  0E0 

i .  6oooo£o 

1  .  8  0  0  0  OE 0 

2  .  0  0  0  0  OE 0 

S2+SA VE 


Y 

1 . 98  669E’“l 
3 . 8941 8E~ 1 
5  .  6  464  3E~ 1 
7 . 17356£"l 
8. 4147 1E~ 1 

9 . 32039£"l 
9. 8  5  4  5  OE  _ 1 
9 .  99574E’~1 
9 . 73848£“l 
9 . 09297F"! 


V 

1 . 96013#"! 
1  .  8  421  2£'-l 
1. 65  067£’”l 
1.39341E  1 
l . osoeoff"! 

7. 24716E“2 
3  .  399  34£_2 
5.8399  1E1  3 
4. 5 44 04£’_ 2 
8.32294F  2 


ItfCY)  I 

2  .  54994£’"9 
1  .  63372£’~8 
3 . 76628£"8 
6 . 21940£"8 
8  .  5  29  2  7 E~ 8 

1 . 02330£”7 
1.09024E  7 
1.01759ff“7 
7 . 78560E_8 
3.5801 IE  8 


IFOOI 

1 . 4982 1E~ 8 
2 . 64235£’"8 

3 . 3  4  7 4 551- 8 
3  .  5  5  99  1E~ 8 
3. 25  873£’~8 

2 . 4  6  01  02?”  8 
1 . 21  349£’_8 
3 . 99  562E’~9 
2.27021F  8 
4. 2  6  8  4 4£’~ 8 


SIGDIG  4 

6 

VYEXACT CD]V 

V  M+D  YEXACT  X 
Cl]  M*-Dp  0 

[2]  M-*-(  SIN  X)  ,COS  X 

V 

RATIO  X 

RATIO  OF  ERRORS 

E(R1)/E(B 2)  E(R1)/E(R 2)  E(B2)/E(R2) 


Y 

V 

Y 

V 

Y 

V 

1050 

2.503 

1045 

1.186 

0.9954 

0.4738 

431 . 7 

17.89 

317  .  5 

9  .  33 

0 . 7354 

0. 5216 

261.  3 

34 . 42 

192.7 

20.06 

0. 7376 

0.5828 

176.6 

53.04 

138.2 

35 . 45 

0 .7825 

0.6684 

122.7 

75.25 

105.1 

60.54 

0. 8564 

0 . 8046 

82  .  75 

103.7 

80.02 

111  .  7 

0.967 

1.077 

49 . 76 

143.9 

56.67 

290.2 

1.139 

2 .017 

19.97 

210.1 

28.72 

1057 

1.438 

5  .  032 

9.18 

353  .9 

19.25 

210.9 

2  .  097 

0.  596 

40.03 

1027 

193.2 

120.6 

4.828 

0. 1174 

4 


SIGDIG  6 


■ 


. 


, 


. 


100 


(5.10) 

VY1PRIME\_U~\V 

V  V+X  Y1PRIME  Y 

[1]  7«-J[2]  ,-Y[l] 

V 

X 

0  0.2  2  4 

Y 

1  0 


USE  1  R 1 
X  5X1  J 


X 

Y 

V 

IS’(X)  | 

\E(V)  | 

2 . 000005“l 

9 . 800675“ 

1 

“3  .  9733351- 

2 

8  .  882645’“ 

8 

5 . 328235“ 

7 

4 . 000005~1 

9  .  2 1 0  6  2# 

1 

“7  .  788265’“ 

2 

1. 23266 E 

6 

1 .037355 

6 

6 . 0  0  0  0  0E~ 1 

8. 25339 E~ 

1 

1 .  129275’“ 

1 

3  .  3  5  7 7  95’“ 

6 

1 . 451535" 

6 

8.  0  0  0  0  05“  1 

6 . 9671 35 

1 

“1.434695’ 

1 

6  .  310325*“ 

6 

1 .718915 

6 

1 . 0000050 

5 . 403125 

1 

1 .68292 E 

1 

9  .  8  6  49  851 

6 

1 .792405“ 

6 

1 . 2  0  0  0  0#  0 

3. 62371 E~ 

1 

“1 .  864065’“ 

1 

1.  373865’“ 

5 

1 . 637645“ 

6 

1 .4000050 

i. 699855“ 

1 

1 .970895’ 

1 

1 .  760675’“ 

5 

1 .235635 

6 

i . 6000050 

"2  .  917845“ 

2 

“1.  9  9  91 45? 

1 

2  .  1  1235 E~ 

5 

5 . 844835“ 

7 

1  .  8  0  0  0  OE 0 

“2 . 27178 E 

1 

“1 .  947705’“ 

1 

2.  394355’“ 

5 

2.997925“ 

7 

2.0000  OE 0 

“4. 16121 E 

1 

“l .  siseis1- 

1 

2.  574275’“ 

5 

1 . 383535“ 

6 

S+SA  VE 

VY2PRIMEIU1V 

V  7«-X  Y2PRIME  Y 

]  V+-YI 1] 

V 

USE 2  B  2 

X  RK2  Y 

X 

Y 

V 

l£(y) 1 

15(7) | 

2 . 0  0  0  0  05“ 1 

9 . 80067 E~ 

1 

"3.97339  S’" 

2 

6. 3  4  6  3  95’“ 

11 

1.  520995“ 

10 

4. 000005  1 

9 . 21061 E 

1 

7. 78837 E 

2 

7 . 797625’“ 

10 

1.705865 

9 

6.  0  0  0  0  0E~ 1 

8  .  253365 

1 

“1.  12  9  2  85’ 

1 

3  .  5  0  5  62E 

9 

4. 545575 

9 

8 . 0  0  0  0  0E~ 1 

6 .967075“ 

1 

"1. 4  3  47  IE 

1 

9. 419425’“ 

9 

8 . 448925 

9 

1 . 0000050 

5 .40302 E~ 

1 

“1.  682945’“ 

1 

1 .950485’ 

8 

1  .  310365 

8 

1 . 2000050 

3  .  623585’“ 

1 

“1.  864085’“ 

1 

3 .430 37E~ 

8 

1 . 812275“ 

8 

1 . 4  0  0  0  OE 0 

1  .  6  9  9  6  75^ 

1 

“1.  970905’ 

1 

5 .401  OSS’- 

8 

2 .306955“ 

8 

1 . 6000050 

2 .  gages'- 

2 

“i.  gggiss1 

1 

7  .  8  3  6  2  451 

8 

2 . 748465 

8 

1 . 80000E0 

"2 . 2720251- 

1 

“1 .  947705”“ 

1 

1.  066375’“ 

7 

3 . 091405“ 

8 

2 . 0000050 

4. 161475’“ 

1 

“i.  siseos1- 

1 

1.376785’“ 

7 

3.293895 

8 

5 1«-S,4  75 


V 

.  . 


* 


v  ,  -  . 

■ 


101 


(5.10) 

USE 2  R2 
X  RK2  Y 

X 

2  .  00000#~1 
4 . 0  0  0  0  0#_ 1 
6.00000#  1 
8  .  0  0  0  0  0E~ 1 
1 . 00000#0 

1 .  20000#0 
1 . 40000#0 
1 . 60000#0 
1 . 80000#0 
2  .  0  0  0  0  OE 0 

S2+-SA  VE 


Y 

9  .  8  0  0  6 7#”  1 
9 . 210  6  IE _1 
8  .  25336#_1 
6  .  96707#_1 

5 . 4  0  3  0 2#~1 

3  .  6  2  3  5  7#~1 

1  .  69967#_1 
2 . 92000#_2 

2  .  2  7  2  0  3E  _1 
4.16147 E  1 


V 

3 .97339E~2 
7 . 7  8  8  3 1E~  2 
1 . 12929E~1 
1 . 43471#"l 
1.6 829 4#  1 

1  .  8  6  4  0  8#~  1 
1 .  97090#~1 
1 .  99915#*! 
1.  9  47  7  OE _ 1 
1.81860#  1 


|#( X) I 

1  .  7  8  39  0#~  8 
4 . 96897#_8 
9.62414#  8 

1  .  5  6  91 6#~  7 

2  .  2  9  9  2  0E~  7 

3. 121 3  5#*7 
3.99401#  7 
4 . 86648#*7 
5. 68169#_7 
6.37944#  7 


l#O0 | 

1 .  43115#*8 
3. 03202#_8 
4.67580#  8 
6.22732#  8 
7 .  5  5  152#*8 

8 . 52143#*8 
9 . 0  26  3  7#~  8 
8  .  9  7  91 8#~  8 
8.32240#  8 
7.0  328  8#” 8 


SIGDIG  4 

6 

7YEXACT [□] V 

V  M+D  YE X ACT  X 

[1]  M+DpO 

[2]  M-*-(  COS  X)  ,  -SIN  X 

V 


RATIO  X 


E(R1 )/E(B2 ) 
Y  V 


RATIO  OF  ERRORS 

#(#1  )/#(#2  ) 

Y  V 


#(B 2  )/#(# 2  ) 
Y  V 


1. 400#3 
1  .  581#3 
9  .  5  7  8#2 
6  .  6  9  9#2 
5 . 0  5  8#2 


3 . 5  0  3#3 
6  .  0  8 1#  2 
3.  1 9  3#2 
2 . 0  3  4#2 
1 . 3  68#2 


4  .  9  7  9#  0 
2 . 481#1 
3. 48  9#1 
4 . 021#1 

4. 2  9 1#1 


3. 7  2  3#1 
3 , 42 1#1 
3 . 1 04#1 
2.  7  6  0#1 
2.  3  7  4# 1 


3  .  5  5  8#_3 

1 .5  69#  2 
3 . 643#"2 
6 . 003#"2 
8. 483#~2 


1 . 063#  2 
5 . 6 2  6#~2 
9 . 721#~2 
1 . 357#_1 
1 . 735#“l 


4 . 0  0  5#2 
3 . 260#2 
2  .  6  9  6#2 

2 . 2  4  5#2 
1 . 870#2 


9 . 0  3  6#1 
5 .  3  5  6#1 
2 . 1 2  7#1 

9 . 6  98#0 

4 . 2  0  0#1 


4 . 4  0 1# 1 
4 . 408#1 

4. 3  4 1#1 
4.21 4# 1 
4 .035#! 


1 . 922#1 
1  .  3  69#1 
6  .  5  0  9#  0 
3 . 602#0 
1  .  96  7#1 


1 . 099#_1 
1 . 352#_1 
1 . 610#_1 
1 . 877#_1 
2 . 15  8#“l 


2. 127#_1 
2 . 556#"1 
3 . 061#”1 
3. 7 1 5#~ 1 

4. 6  8  4#~ 1 


4 


SIGDIG  6 


. 


■ 


